Hall viscosity, orbital spin, and geometry: paired superfluids 
and quantum Hall systems 
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The Hall viscosity, a non-dissipative transport coefficient analogous to Hall conductivity, is con- 
sidered for quantum fluids in gapped or topological phases. The relation of the Hall viscosity to the 
mean orbital spin per particle s (discovered in previous work) is elucidated with the help of exam- 
ples and of the geometry of shear transformations and rotations. For non-interacting particles in a 
magnetic field, there are several ways to derive the result (even at non-zero temperature) , including 
standard linear response theory. Arguments for the quantization, and the robustness of s to small 
changes in the Hamiltonian that preserve rotational invariance, are given. Numerical calculations of 
adiabatic transport are performed to check the predictions for quantum Hall systems, with excellent 
agreement for trial states. The coefficient of k 4 in the static structure factor is also considered, and 
shown to be exactly related to the orbital spin and robust to perturbations in rotation invariant 
systems also. 
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I. INTRODUCTION 

In recent work the notion of Hall viscosity in quantum 
fluids [l| has been revived @, Q , and values of this pa- 
rameter have been calculated for several systems p| |1] 
(the earlier works termed it "odd" or "antisymmetric" 
viscosity [l|, or "Lorentz shear modulus" Q; the term 
"Hall viscosity" is from Ref. 3]). The underlying def- 
inition can be briefly described as follows [l(: For an 
elastic solid in d dimensions (d = 2, 3, ...), the low- 
energy, long- wavelength effective stress tensor <7 b(x, t) 
determines the local force density f a on the system by 
f a = — J2b d(Jab/dxb. The stress tensor can be expanded 
in powers of the local strain u a &(x, t) (relative to a relaxed 
or unstrained configuration) and its derivatives Q: 
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where X a bef is the tensor of elastic coefficients (moduli), 
and rjabef is the viscosity tensor (a, 6, . . . = 1, . . . , d). 
In an isotropic solid, a and u are symmetric, so we have 
Xabef — Xbaef = X a bfe, and the same for r\. The linearized 
strain is given in terms of the displacement u a (x, t) from 
the unstrained configuration by 



U ab =\ 



dUg 

dx b 



du b 
dx a 



(1.2) 



Similarly, for an isotropic fluid, the displacement from the 
arbitrary choice of unstrained configuration should not 
enter, so the elastic moduli vanish, with the exception of 
a bulk term in a a b which is p5 a b where p is the pressure (p 
depends on the density). In place of the time derivative 
of the strain, one has the symmetrized derivatives of the 
velocity field v(x, t), 
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In addition, for a fluid the momentum flux /iv a Vb (where 
/i(x, t) is the mass density) must be included as part of 
a ab ; then dg a /dt + J2b d(J ba/dx b = 0, where g(x,t) is the 
momentum density. 

Now (in either a solid or a fluid) r\ can be divided 
into symmetric and antisymmetric parts with respect to 
interchanging the first with the second pair of indices; 
thus 7 labef = r) { a fl f + rf^l } , where 
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Only the symmetric part contributes to dissipation of en- 
ergy, as can be seen from the rate of entropy production 
(per unit volume) ds (x, t)/dt due to the above stress ten- 
sor [1, 
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(where j s (x, t) is the entropy flux), and r/^ should be a 
positive quadratic form. In a gapped quantum fluid at 
zero temperature, this part should vanish. The antisym- 
metric part t/' 4 -' is termed here the Hall viscosity tensor. 
It is a set of non-dissipative (or reactive) transport co- 
efficients, and is closely analogous to the antisymmetric 
Hall conductivity tensor. The analogy is best seen by 
viewing the Hall viscosity as the stress response to an 
applied field, which here is a time-dependent metric ten- 
sor, as we will discuss shortly. In view of this analogy, 
there should be no more objection (on the grounds that 
it is non-dissipative) to terming the former a viscosity 
than there is (on the same grounds) to terming the latter 
a conductivity. It is not an elastic modulus as the corre- 
sponding stress vanishes in the static limit, and is not in 
general connected with the Lorentz force. 

In an isotropic fluid in d dimensions, the symmetric, 
dissipative viscosity tensor is determined by only two co- 
efficients, the familiar bulk and shear viscosities (which 
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must be non- negative). For the Hall viscosity, rotational 
invariance forces rf A > to vanish identically if d > 2, but 
in two dimensions there is a single rotationally invari- 
ant tensor, and we denote the corresponding coefficient 
by rjW also Q. On the other hand, even a mildly non- 
isotropic fluid in more than two dimensions (for example, 
one in which rotational symmetry is spontaneously bro- 
ken by the appearance of an intrinsic angular momentum 
only) can have a non-zero Hall viscosity, and we will have 
more to say about that later in this paper. However, we 
note also that the Hall viscosity is odd under both time 
reversal [l| and (in two dimensions) reflection of space, 
and so must vanish when either of those symmetries is un- 
broken (or in the case of reflections in non-rotationally in- 
variant systems, some components corresponding to un- 
broken reflection symmetries must vanish) . Hall viscosity 
has been known for some time in classical plasmas in a 
magnetic field @. 

Avron, Seiler, and Zograf (ASZ) [1] related the Hall 
viscosity to the adiabatic response to a slowly varying 
metric tensor. This parallels certain formulations of Hall 
conductivity as a Chern number or as adiabatic response 
and in the present case is based on the fundamen- 
tal fact that varying a Hamiltonian (or an action) with 
respect to the metric tensor produces the stress tensor. 
ASZ calculated the Hall viscosity of the filled lowest Lan- 
dau level (LL) in the non-interacting case, and found that 
it is an intensive quantity, independent of the shape of 
the fluid. Independently at around the same time, Levay 
calculated the same adiabatic curvature for a single par- 
ticle in any LL jl| . This may be used to extend the ASZ 
result to more general filling factors, and also (by per- 
forming a thermal average) to recover the classical result 
at high temperature [3J. 

In the recent work by one of the authors [H , the adia- 
batic approach to Hall viscosity was generalized to some 
other systems, mainly in two dimensions: paired (and 
gapped) superfluids, and fractional quantum Hall (QH) 
wavefunctions, starting with the Laughlin states [l(j. It 
was realized that the Hall viscosity can be written in the 
form 
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where n is the particle number density, and s can al- 
ways be naturally interpreted as minus the mean orbital 
spin per particle (this is not always the total angular mo- 
mentum per particle). For non-interacting particles in a 
magnetic field, this spin is due to the cyclotron motion. 
For paired states, it is the intrinsic angular momentum 
of a Cooper pair. This s is also related to the "shift" S, 
an offset that is required in the number of magnetic flux 
quanta (in units of hc/e) piercing the surface when the 
ground state is formulated on a sphere: 



S = 2s. 



(1.8) 



Both properties are expected to be robust (quantized) in 
a quantum fluid as long as translation and rotation in- 
variance are not broken. We note that the notion of an 



orbital spin was invoked in Ref. 11] in order to explain 
the shift. Two subsequent pap ers have rederived the re- 
sult for the Laughlin states [12| , and attempted a general 
discussion [13l |. 

In the present paper, our main goals are to add some 
insight into the general picture just described, and to 
present numerical tests of the results of Ref. Q . In par- 
ticular, we wish to explain why the stress in the adiabatic 
response to a strain is related to some sort of angular mo- 
mentum. The general theory of adiabatic response, which 
we review in the main text in Sec. [HI states that the stress 
can be obtained by varying the Hamiltonian with respect 
to the metric, or by applying a strain. A second variation 
gives the stress response to a slowly varying strain. Now 
a small uniform strain (or uniform change in metric) that 
preserves the area can be described by a symmetric ma- 
trix with constant real coefficients, which is close to the 
identity and has determinant 1; typical examples in two 
dimensions are 
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where e, e' are small. The first of these stretches the 
x coordinate, and squashes the y coordinate. The sec- 
ond does the same but for axes rotated by 7r/4. Then 
the adiabatic response theory relates the response to the 
commutator of the effects of two such changes (which de- 
scribes a Berry phase). The effects of applying two such 
strains in opposite orders differ by a small rotation: 
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If the system is in an angular momentum eigenstate, then 
the effect of this rotation is to multiply it by a phase. 
Thus the Berry phase is related to the spin of the system. 
We will go through this argument in much more detail in 
what follows. In practice, we need to consider the system 
on the torus (i.e. periodic boundary conditions), and it 
may not be an exact eigenstate of angular momentum, 
but we will show that the above simple picture captures 
the essence of the situation. As far as possible, we use an 
operator approach, and a type of vector bundle called a 
"homogeneous bundle" (related to spin coherent states), 
that generalizes the approach of Ref. [9( . 

In addition, we present in Sec. [IT] various relevant exam- 
ples of calculations in toy models, including trial states 
in the fractional QH effect in the disk geometry. We 
give the explicit generalization to more than two dimen- 
sions. We show that the result for non-interacting par- 
ticles in a magnetic field can be reproduced by standard 
(Kubo formula) linear response theory. We present ar- 
guments for the quantization of s to values that are ra- 
tional numbers, and for its robustness under perturba- 
tions of the Hamiltonian subject to the requirement of 
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translation and rotation invariance (apart from bound- 
ary conditions). We show numerically that the predic- 
tions of Ref. [3] for the Hall viscosity of some fractional 
QH trial states [the Laughlin and Moore- Read (MR) [l4| 
states] are correct. We show that the results are robust 
to a perturbation of the Hamiltonian in one case, and 
examine the effect of a particle-hole symmetry-breaking 
three-body term on the result at half filling of a LL. We 
propose that the adiabatic calculation of s can be used to 
determine the shift in the ground state at a given filling 
factor, and thus to obtain information about the topo- 
logical phase in question. 

In a separate approach, motivated by some remarks in 
Refs. [Hill, we consider in Sec. Hill the static structure 
factor of a QH ground state, and the related compress- 
ibility xt of an equivalent two-dimensional (2D) system. 
We find exact results for \t and for the coefficient of 
the k term in the structure factor, under a hypothesis 
similar to one used in Ref. [3j . This can be tested numer- 
ically by Monte Carlo simulation, with good agreement; 
the full static structure factor s(k) for the MR state is 
also obtained. We find analogous results also for wave- 
functions for anyons in zero magnetic field, where we find 
Xt and s(0). For QH states in the lowest LL, the exact 
Xt allows us to recover the Hall viscosity again (T^ |. 

II. HALL VISCOSITY 



bundle, which will be useful in the calculations. In Sec. 
Ill A 5[ we give an explicit example of such a homogeneous 
bundle, based in single-particle quantum mechanics, and 
show how the adiabatic curvature is obtained. In Sec. 
Ill A 6[ we introduce periodic boundary conditions, and 
point out that while the resulting bundle is not a homo- 
geneous bundle, it can be well approximated as such in 
some limits. With this we then recover the Hall viscosity 
of paired superfluids in two dimensions, which is related 
to the angular momentum of the pairs; the result is the 
same as in Ref. [3|, but the derivation differs. In Sec. 
Ill A 71 we turn to 2D systems in a magnetic field, starting 
with a single-particle model as before. We discuss differ- 
ent ways to deform the system, before recovering the Hall 
viscosity of non-interacting particles as in 1,9]. We show 
that the result can also be obtained by direct Kubo linear 
response theory. For correlated or fractional QH states, 
the result obtained in Ref. [3[ is recovered, but now using 
the disk geometry. In Sec. Ill A 81 we briefly discuss the 
situation for paired superfluids in three dimensions. In 
Sec. Ill A 91 we discuss the issue of the robustness of the 
spin per particle s, to perturbations of the Hamiltonian. 
We first give one very direct and compelling argument. 
Then we also discuss a brute force approach in perturba- 
tion theory. Finally, we also discuss robustness in terms 
of the claim that 2s is the shift, which is a rational num- 
ber. Some further side discussion related to this Section 
appears in Appendices [B] and [Cj 



In this Section, we discuss the Hall viscosity. Sec. Ill Al 
discusses all the analytical arguments in some detail. Sec. 
Ill B I discusses our numerical work, in which the analytical 
predictions are tested. 



A. Analytical approach 

We will give a direct a priori derivation of the rela- 
tion of the adiabatic curvature to the total orbital spin, 
and hence of the Hall viscosity with the orbital spin den- 
sity, noticed in [3]. We will be satisfied with the case 
of a system with translational and rotational symmetry, 
for which the quantities are quantized, though the ap- 
proach is more general. However, we do discuss higher- 
dimensional situations that are not truly isotropic. 

In more detail, in Sec. Ill A II we review the general 
set-up for calculating linear response of a system to 
a perturbation to adiabatic transport, using the no- 
tions of adiabatic or Berry connection and curvature 
(Berry phase) . This approach is generally useful for non- 
dissipative transport coefficients in a gapped system. In 
Sec. Ill A 21 we introduce notation for some coordinate 
systems, and for describing deformations of the metric. 
We also define some group-theoretic notions. In Sec. 
Ill A 31 we combine the preceding ideas to arrive at ex- 
plicit expressions for the Hall viscosity in terms of adi- 
abatic curvature (following ASZ [1]). In Sec. Ill A4l we 
introduce the general idea of what we call a homogeneous 



Adiabatic response and transport 



This Subsection, Sec. Ill A II is a short review of adia- 
batic transport [l5l - [l7j and response theory p| [l8j]; the 
results will be applied to Hall viscosity afterwards. We 
suppose that the Hamiltonian H (A) depends on a set of 
parameters ( "generalized coordinates" ) A = {A M } (// = 1, 
. . . , n), and we are interested in a particular eigenstate 
|<p(A)) of H(X), which we assume for the present is non- 
degenerate, and also is separated by a gap from the rest 
of the spectrum of H(X); these statements should hold at 
least in a neighborhood of the origin A = (thus no lev- 
els cross). By subtracting the energy eigenvalue of \(fi(X)} 
from -ff(A), we can assume that this eigenvalue is zero, 
for all A. 

Our interest is in calculating the linear response of 
some "current" or "generalized force" in the state \ip) 
to the application of some "field" or "generalized veloc- 
ity". We assume that current operators of interest can 
be written as 



(2.1) 



We will also write = d/dX^. Then the current in a 
state \if}) is 



I M (X) = W M (A)|</>). 



(2.2) 
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For = \<p(X)), it follows from H(X)\(f(X)) = for 
all A that / M (A) = 0, which expresses the absence of 
"persistent currents". The field or generalized velocity 
in the generalized coordinates is represented by a time- 
dependent A, with generalized velocity V„ = dX v /dt. A 
basic example in the QH effect is that in which the sys- 
tem has periodic boundary conditions (i.e. is topologi- 
cally a torus), and the coordinates A are the Aharonov- 
Bohm fluxes (line integrals of the vector potential) <j>i, 
</> 2 through the two cycles of the torus. Then 7 M is the 
current, and V„ is minus the electric field. A weak field 
corresponds to slow variation of A, and so linear response 
will correspond to adiabatic response to time-dependent 
A. In the present case, the fact that we wish to take a 
derivative with respect to A^ means that we must con- 
sider adiabatic transport at any A in some neighborhood 
of the origin, and so we use a velocity field V„(A), which 
depends on A but not explicitly on t. 

The quantum adiabatic theorem asserts that if 
\ipx>(X)) is the normalized solution to the time-dependent 
Schrodinger equation 

ff(A)|Vv(A)> = ift|Vv(A)), (2.3) 

(where A = \(t) along an integral curve of V„ is under- 
stood, and note that dt = ^2, v Vvdv) with initial con- 
ditions \ip) = \lp{\)) for A on some n — 1-dimensional 
surface transverse to V„ at time t — 0, then the limit 
IV'vPO) -> |V>P0) ( as Vv -> 0) exists, and \ip(X)) is a 
multiple of \ f{X)} (in general, \ip(X)} is not independent 
of V V (X) from which it was obtained). Hence (by consid- 
ering different initial values of A on the surface) we have 
a smoothly varying state \il>(X)) 10 r all A. 

It follows from the adiabatic theorem [l5|, [T(| that in 
the limit, the state \i/j(t)) at any A obeys 

(m^) = 0, (2.4) 

which means that is parallel transported by the adi- 
abatic evolution along the curves of V v . Formally, we 
can view this set-up as a vector bundle over the manifold 
with coordinates A, in which the fibre at each point A 
is the one-dimensional vector space spanned by |<p(A)). 
In terms of the basis states \ip(X)), there is a Berry or 
adiabatic connection (vector potential in A space) on the 
bundle, given by 

A„=i{i P \d ll <p). (2.5) 

For the adiabatic linear response, we differentiate eq. 
(j2.3[) with respect to A M , to obtain 

d^H\ih) +H\d^ v ) =iJ2 d M d M + id t \dM- 

(2.6) 

Using this and eq. (12.31) itself in eq. (I2.2[) . we obtain 
J M (A) - -idt^vld^v) -iJ2 d M4>v\d^v) (2.7) 



[which appears to differ from eq. (5.12) in Ref. [IH], and 
then finally taking of eq. (|2.4p for adiabatic transport, 
we obtain, to linear order in V v , 

ip(A) = iJ2 [9M d vi>) - dA^M V„. (2.8) 

The expression in closed brackets is gauge-invariant 
(i.e. it is invariant under multiplication of \tp) by a A- 
dependent phase factor), and so the arbitrary smooth 
basis \ip) can be used instead of the parallel-transported 
| ip) . Then we have 

i P (A) = 53F |w (A)V„(A), (2.9) 

V 

where 

= d„A v - dvAp (2.10) 

is the adiabatic curvature or field strength of the connec- 
tion A above; it is manifestly gauge invariant, as claimed, 
and hence independent of V„. This final expression for 
the non-dissipative transport coefficients, given by F, can 
also be obtained by other approaches based on conven- 
tional Kubo linear response theory, rather than adiabatic 
transport, as has been shown explicitly for the Hall con- 
ductance 

More formally, the set-up for a g eneral Berry or adia- 
batic connection calculation [la - flTf involves again a vec- 
tor bundle embedded in a Hilbert space %. A bundle has 
a base space that is a manifold, again with coordinates 
A = {A M }. Over each point A, there is a fibre V\, which is 
a subspace of T-L, with the same dimension (referred to as 
the dimension of the bundle; it can be finite or infinite) 
at all A; we write V for V\ at generic A. A continuous 
vector function \ f{X)) of A with |y(A)) G V\ for all A is 
called a section of the bundle. There is an inner product 
on each fibre V\, determined by that on H. We take or- 
thonormal basis vectors \tp a {X)) (labeled by a — 1, ... , 
dim V) in each V\ , that are a collection of sections which 
vary smoothly with the coordinates A. Then the Berry 
or metric connection [l6[ is a Hermitian matrix function 
of A given by [H-Qj]] 

A^ap = i{<Pa\dn<Pp)- (2-11) 

In components, a section is \<f(X)) = J2 a v a (X)\ip a (X)) , 
and the covariant derivative is then \D,j,<p(X)) = 
E a/ ?(^^^ - iA^. a pv^)\ip a (X)). Parallel transport of 
a vector along a curve with tangent is defined by 
the condition £ M V/JD^A)) = 0, or by ^ V^d^Va - 
i^2p A^_ a pvp) = in components. Parallel transport 
around a small loop picks up the integral of the curva- 
ture of the connection, given in matrix notation by 

Ffj, v = i[D„,D u ] (2.12) 
= d^A v -d u A^i[A^A v ]. (2.13) 

For a one-dimensional fibre the curvature reduces to 

F^y = i {{d^dyip) - {d v ip\d^)\ (2-14) 
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as above. This curvature or field strength for dimT^ > 
1 is covariant under a unitary gauge change (varying 
smoothly with A) of the choice of orthonormal basis for 
V\. For the applications here, the fibre may be spanned 
by a single vector, or may have dim^ > f. (A vec- 
tor bundle in which the fibre is one-dimensional is also 
called a line bundle.) We emphasize that the connec- 
tion (covariant derivative) and curvature depend on the 
choice of the fibre subspaces V\ of the Hilbert space. As 
an extreme case, if V\ is the whole Hilbert space rl for 
each A, then an easy calculation from eq. (|2.f 3j) (using 
\ faW)(faW\ = I) shows that F^afs = 0. This is 
because \cp a (X)) differs by a (unitary) gauge transforma- 
tion from a A- independent basis for H, so F^ap must 
vanish. 



2. Deformations of shape and of metric 

We now consider some geometric aspects of the prob- 
lem of uniform deformations of an object. A general 
linear transformation of the Cartesian coordinates in d- 
dimensional space (d > 1) that leaves the origin fixed is 
described by an invertible d x d matrix with real entries, 
and so is an element of the group GL(d, R). Thus writing 
a typical element as a matrix A, it acts as 



( X1 \ 






X2 








— > A T 











onx= [x\, . . . , Xd) and this is viewed as an active trans- 
formation. (The transpose T of A is used to make later 
notation simpler.) We will consider only the transforma- 
tions that can be continuously deformed to the identity 
(and not e.g. reflections); these matrices have positive 
determinant, and form the connected group GL + (d, R). 
For all d, this group contains a subgroup that consists 
of the positive multiples of the identity, and is isomor- 
phic to the group R+ of positive real numbers (under 
multiplication) ; these elements represent the simple scale 
transformations or dilatations. As we intend to study 
mainly incompressible systems, we will usually restrict 
to transformations that preserve the volume, which are 
represented by matrices A with det A = 1; the group of 
these is denoted SL(d, R) (we also restrict to d > 2 from 
here on). In any case, even if the subgroup = R* is 
included, it decouples from the following considerations 
because its elements commute with all others. 

We now fix a coordinate system with coordinates x a 
(a = I, . . . , d). This will be important in setting up 
quantum mechanics, which will work in a fixed Hilbert 
space, with the norm-square of a state given by integrat- 
ing the square of the absolute-value of the wavefunction 
in these coordinates. We will often write 

X = A T x (2.16) 



when A is given. We assume that the metric for the 
(i-dimensional space has the canonical form in terms of 
displacements of X. That is, ds 2 = ^2 ab S a bdX a dXb = 
J2ab9abdx a dx b , here always with g independent of posi- 
tion. Then the metric, viewed as a matrix, is 

g = AA T . (2.17) 

Linear transformations change the metric by g — > A'gA' T , 
but leave the x coordinates unchanged. Thus we have 
parametrized metrics by transformations A € SL(d, R) 
from an initial metric which we take to be g = I (7 = I n 
is the n x n identity matrix). 

The group G = SL(d, R) contains a family of compact 
subgroups isomorphic to SO(d). In particular, there is 
the subgroup (which we will denote by K) consisting of 
d x d orthogonal matrices O € G obeying O t = O -1 . 
These leave the initial metric g = I invariant, and so 
any g is invariant under A — > AO for any O G K. Con- 
sequently, the possible metrics g are in one-to-one cor- 
respondence with the cosets AK, which are the points 
of the coset space G/K = SL(d,R)/SO(d). The latter 
is a classical example of a Riemannian symmetric space 
[l9l [20| . An arbitrary metric g — AA T is invariant un- 
der g — !• OgO T for O of the form O = AOA^ 1 for all 
O G K. For each A, the set of such O is a group K\ that 
is isomorphic to if. As A varies, these form a family of 
compact subgroups in G. 

Now we pass to the Lie algebra. The Lie algebra ele- 
ments a correspond to elements A £ G by A = e a . Then 
the Lie algebra of G is represented by the real dx d ma- 
trices with zero trace, and is denoted sl(d, R). The Lie 
algebra of K is represented by the antisymmetric matri- 
ces, and is denoted so(d); it is a Lie subalgebra. Hence 
the quotient space of these two Lie algebras can be viewed 
as consisting of the real symmetric matrices. The latter 
do not form a Lie algebra [because so(d) is not a Lie 
ideal in sl(d, R)]. Instead, the commutator of two sym- 
metric matrices is an antisymmetric matrix (an example 
was shown in Sec. [TJ) . This means that if we consider two 
infinitesimal shear transformations of the metric g = I, 
represented by symmetric matrices, their commutator is 
antisymmetric and generates a rotation. Further, the 
commutator of an antisymmetric with a symmetric ma- 
trix is symmetric. [There is a similar picture for any g, 
but in terms of K\ and the quotient of sl(d, R) by the 
Lie algebra of K\.) The decomposition of the Lie alge- 
bra into two parts in this way, in which AT is a maximal 
compact subgroup, is called the Cartan decomposition, 
and is the structure arising in all Riemannian symmetric 
spaces [l9L [20([ . 

We add here some additional background on the struc- 
ture of the Lie groups. SL(2,R) is closely related to 
SO(2, 1), the group of real linear maps preserving a sym- 
metric bilinear form of signature (2, 1). As in the better- 
known related case of the corresponding compact groups 
SU(2) and SO (3), SL(2,R) is a double cover of SO(2, 1) 
(there is a two-to-one group homomorphism from the for- 
mer to the latter, with kernel Z2 = {±72}). However, 
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unlike the compact versions, SL(2,R) is multiply con- 
nected [like SO(2)], with fundamental group 7r 1 (G) = Z 
(and hence 7Ti = Z for SO(2, 1) also). Consequently, it 
admits M-fold covers for all integers M > 0. There is a 
universal cover of SL(2, R) that is simply connected, but 
it cannot be faithfully represented as a group of finite 
matrices. For d > 3, SL(d, R) is doubly connected [like 
SO(d) for d > 3], so it has a simply-connected double 
cover. For all d > 2, we will denote the double cover of 
SL(d, R) by SL(d, R). We emphasize that finite covers 
of a Lie group always have the same Lie algebra; locally, 
they are the same groups. 



In terms of A, 

5 9ab = [(^A.A- 1 ).^ + g ac (5AA- l ) bc ] . (2.20) 

c 

Now SA.A^ 1 is rather complicated in terms of the global 
coordinates A, A = e x . It will be simpler to use local 
coordinates A' in the vicinity of any given g, defined by 
the left action of G on A (sometimes called left translation 
on G), that is A -> A' A for A' = e A ' e G, with A' small 
(A' = at the given g or A). Then SA.A^ 1 = 8X', and 



lab 



= J2( SX 'ac9cb + gacS\' bc ). (2.21) 



3. Expressions for Hall viscosity 



The parametrization in Sec. Ill A 21 is useful for deriv- 
ing the Hall viscosity, because it provides a system of 
non-redundant coordinates A, and the use of such coor- 
dinates was assumed in the formulas of Sec. Ill A II (The 
condition that the matrices A be traceless is a remain- 
ing constraint, but we will eventually see that it can be 
dropped.) For the continuity equation involving momen- 
tum and stress of a many-particle system to hold, the 
Hamiltonian H must be translationally invariant. To ob- 
tain the stress tensor by adiabatic response, we begin by 
noting that it can be defined as 



0~ab 



s— , 

5g a b 



(2.18) 



where again H is the Hamiltonian, and S/Sg is a func- 
tional derivative, applied by varying g at a position x 
only. (This may be checked by considering the mo- 
mentum flux part, which arise from the kinetic energy.) 
This definition produces a symmetric tensor, because 
the metric tensor is symmetric. The symmetry of the 
stress tensor follows from rotational invariance. In a non- 
rotationally-invariant system, the stress tensor does not 
have to be symmetric, and is not given by this formula; 
the symmetry analysis of Section U does not hold in this 
case. 

To obtain the Hall viscosity, we wish to consider a sys- 
tem of particles with periodic boundary conditions on a 
square or cube in x space, and make a uniform variation 
of the metric [![. For a uniform variation of g, we can 
instead write 



Eab — 2 



dH 

dg a b 



(2.19) 



where (throughout the discussion) g is constant in x 
space. Thus £ = L d a, where throughout this paper L d 
will be the volume of a finite system in dimension d. For 
such an expression, in which the symmetric tensor g is 
varied to obtain the symmetric E, we should be care- 
ful because the components of g are not all independent. 
The correct result is obtained by viewing the variation of 
g as unconstrained, and then symmetrizing the tensor of 
partial derivatives. 



Assuming A' enters the Hamiltonian only through g, we 
have 



E q6 = g" 



dH 
dX> 



(2.22) 



where g ab is the inverse metric of g ab , so g g b . 



(Other than for g ab , we will not distinguish upper and 
lower indices, but we may note that if we did, then the 
coordinates x a , the stress a ab as defined above, and vis- 
cosity r\ a bcd would each have all indices upstairs, while 
momentum density g a has a down index, and A has first 
index down, second index up. Indices can be raised us- 
ing g ab or lowered using g ab . The continuity equation for 
momentum is then dg a /dt + da a b /dx b — 0. The distinc- 
tion of up and down is ultimately unimportant as we set 
g = I at the end of the calculation.) The right-hand side 
of eq. (|2.22p will be automatically symmetric, but can be 
explicitly symmetrized if there is any concern over this. 

Similarly, for the definition of viscosity, we can identify 
duef/dt = ^dg e f/dt, and write 



du 



bef 



ef 



ef 



dt 



/ d\' d\' 



efc 



^2 VabefQe 
efc 



V dt 

dt ' 



(2.24) 



using the symmetry of n a bef under ef>/. 

We may now compare these definitions with that for 
adiabatic response, in which the indices like /i will be 
replaced with pairs like ab, which appear on the non- 
redundant local coordinates X' ab . We find (assuming a 
non-degenerate state for the present) 



T]abcd 



be 



9 df F ae , cf . 



(2.25) 



ef 



This tensor r\ — rp J is explicitly antisymmetric under 
exchange of pairs ab O cd, and also symmetric under 
a O b and under c O d if the calculation has been set up 
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as we defined it. If there is any doubt (and the system is 
translation and rotation invariant), symmetrization un- 
der both exchanges can be applied to the right-hand side. 
For a non-rotationally invariant Hamiltonian, A' does not 
enter only through the metric, and moreover the stress 
tensor need not be symmetric. We may nonetheless con- 
sider adiabatic deformations of the shape with A', and 
the result (I2.25|) still applies (if there is a metric), though 
the symmetries of rj discussed in Sec. U are lost, except for 
Vab,cd = —Vr A.ah- In this way we can make contact with 
Haldane [13J. In what follows, we usually concentrate on 
rotationally-invariant systems. 



4- Homogeneous bundles 

Some vector bundles of a particular type will be rele- 
vant here; in mathematics these are called homogeneous 
bundles [21| . Suppose that we have a unitary represen- 
tation W of G in some complex Hilbert space (it need 
not be the entire Hilbert space). As G is non-compact, 
such a representation is either trivial (all elements of G 
act as the identity) or infinite dimensional. We obtain a 
finite-dimensional vector bundle over G/K by first tak- 
ing a finite-dimensional unitary representation Vq of AT, 
with Vb C W (these exist because K is compact). We as- 
sociate this vector space Vq with the origin A = 0, where 
the coordinates A are obtained from any representatives 
for the Lie algebra coset space corresponding to G/K 
[thus A = corresponds to A = / (mod K)]. Then we 
associate a similar vector space with every point of G/K 
by applying the action (in W) of a corresponding A e G 
to Vb- This is well-defined on G/K because Vq is a rep- 
resentation of K. Then we have a vector bundle over the 
base space G/K, in which the fibre V\ over each point 
A = e A (mod K) is isomorphic to Vb- Similarly, by re- 
placing K with the trivial subgroup, we will also consider 
vector bundles over G constructed in a similar way, and 
call these homogeneous bundles also. In all cases, the 
union of the subspaces V\ forms a representation of G, 
and without loss of generality we assume from here on 
that this subspace of Hilbert space is all of the represen- 
tation W. Then W is said to be the representation of G 
generated from the representation Vq of K C G. In many 
cases occurring in our discussion, the representation W is 
irreducible, and its structure can be determined, though 
it is not clear if this information is useful physically. De- 
tailed examples of homogeneous bundles will appear in 
the following discussion. 

A homogeneous bundle is a more general version of 
some constructions that are fairly well-known in physics, 
such as coherent states for a quantum spin (for a review, 
see Ref. [22j and papers reprinted therein); we pause to 
describe these. For coherent states of a compact or non- 
compact Lie group G, W is taken to be a highest-weight 
representation, and the representation Vq is spanned by 
the highest weight vector of W, that is annihilated by all 
the "raising operators" of the Lie algebra of G (the pos- 



itive roots). Then there is a compact subgroup K that 
contains the Cartan subgroup, and which maps the high- 
est weight vector to a scalar multiple of itself. (There is 
a similar construction for lowest weights, of course.) The 
best-known example is G = SU(2), and K the U(l) sub- 
group generated by (say) S z , so G/K is the two-sphere 
S 2 . Then when W is the spin-5 1 representation (S = 0, 
1/2, 1, ...), Vq is spanned by the vector of maximum 
S z = S. Similar bundles, though they are not usually 
called coherent states, can be constructed from the same 
G, AT, and W by choosing Vb to be spanned by an eigen- 
vector of S z with eigenvalue m < S (for m = —S, one 
is back to coherent states constructed from the lowest 
weight). Analogs of these with dimVb = 1 can be found 
for other Lie groups also, by lowering the highest- weight 
vector in a highest weight representation W, with K de- 
fined as the isotropy subgroup that maps Vb into itself. 
In such cases, the representation W is irreducible. (Some 
authors would define coherent states for a group G as any 
example in which Vb is one-dimensional [23[ ; in this case 
W may not be irreducible.) But these coherent-state 
examples are not the most general ones for the construc- 
tion described above. Examples of homogeneous bundles 
with dim Vq > 1 (as we will obtain later) can arise from 
coherent state bundles only as direct sums. 

5. Example of a homogeneous bundle 

Now we will give an example of the set-up described 
in the previous subsections, and demonstrate the basic 
result for the adiabatic curvature (or holonomy, or an- 
holonomy) of a homogeneous bundle over G for G = 
SL(d, R), using this example. For simplicity, we assume 
Vq is spanned by a single vector, which we represent by 
a fixed function /(x) in the Hilbert space H = L 2 (R d ) 
of square-integrable functions, which represents a single 
particle in zero magnetic field. Then V\ is spanned by <p\ 
(we make a small change in notation; this corresponds to 
|<p(A)) in the previous sections), 

p A (x) = /(A T x), (2.26) 

and if a is normalized, / d d x |</?a(x)| 2 = / d d x |/(x)| 2 = 1 
for all A as dct A = 1 . This form could arise for example 
from the family of Hamiltonians 

tf A =--Lv x + *7(X), (2.27) 

ZlTlp 

(where Vx has components d/dX a , m p is the particle 
mass, and U is a potential function) which are to be 
viewed as operators on the Hilbert space of functions 
of x, not X. If /(x) is a non-degenerate eigenstate of 
Hj, then /(A T x) is an eigenstate of H\ with the same 
eigenvalue. Then a short calculation shows that the left 
action of G, that is A — > A'A for A' = e x on <^a as before, 
is given by 

9Kb dx b 
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where the derivatives are taken at A' = 0, and on the left 
hand side x is held constant, while on the right hand side 
/(A T x) is viewed as a function of x with A held constant. 
Therefore we define 

Jab = w«^-, (2.29) 

as operators on %. Their commutation relations are 
those of gl(d, R) (here with factors i included in the gen- 
erators), 

[Jab, Jed] = i(8bcJad — 8 a dJcb)- (2.30) 

The operators J a b are not all self-adjoint, though they are 
for the traceless combinations used to generate SL(d, R). 
It is convenient for us to use A's that are not constrained 
to be traceless, so instead we define 

~ d 

Jab = ix a-Q^T + p^ab, ( 2 -31) 

which are self-adjoint: j\ b = J ab . The traceless combina- 
tions are unaffected by this change — only the generator 
of the subgroup isomorphic to R* of scale transforma- 
tions is affected. The commutation relations of J ab are 
the same as those of J ab . 
Then in W we can write 

|<M)=^ trATJ M. (2-32) 

We define S{A) = e ~ ttlxTj to be the "strain operator" 
that implements the transformation from A = / to A. 
This completes the construction of a homogeneous bun- 
dle over G. In fact, if A is not required to be traceless, 
then the same construction yields a homogeneous bundle 
over A e GL(d, R), in which the wavefunction of <^a(x) 
is 

^ A (x) = (detA) 1 / 2 /(A T x), (2.33) 

and the determinant factor maintains the normalization 
as A varies, as it must because the representation is uni- 
tary (the generators are self-adjoint). However, if we 
want this ip\ to be an eigenstate of H\ above for all 
A e GL(d, R), then U(X.) must in general also depend 
directly on detA, unless U(X) scales as degree —2 (i.e. 
the same as the kinetic term). 

To obtain a homogeneous bundle over G/K in the d = 
2 case [or the same with GL(2,R) in place of G], Vq is 
supposed to be a representation of K, so we assume that 
/ has the form /(x) ~ |/(x)|e~ ls< ^, where x±+ix 2 = re 1 ^ , 
and r > and </> are real (s is an integer). This can arise 
from the family of Hamiltonians H\ if in addition £/(x) 
is rotationally invariant. For the subgroup K, we have 
— J12 + J21 = —id /dip, and \ipi) is an eigenvector with 
eigenvalue — s. We note that for A £ K, A T = A -1 , and 
the definition of <p\ for this coincides with the standard 
action of an active rotation of the state ipi. Hence we 



have obtained a line bundle over G/K for any such choice 
of/. 

The representation W of G associated with such a bun- 
dle has one of two forms. First we note that while the 
eigenvalues of —id /dp are integers, the remaining gener- 
ators of sl(2, R) raise or lower this eigenvalue in steps of 
2, due to the "quadrupolar" nature of the strain. Then 
it turns out that the unitary representation W is irre- 
ducible, and contains eigenfunctions for —id/dpi of all 
possible eigenvalues that equal s (mod 2), each occur- 
ring with multiplicity 1. These representations lie in the 
set known as the principal series of irreducible unitary 
representations of SL(2,R), for which the possible val- 
ues of the quadratic Casimir of SL(2,R) are continuous 
and bounded below (or above, depending on a sign con- 
vention) |24j . In our representations, the value of the 
Casimir depends on /, but will not be needed. 

When we apply the general result for adiabatic cur- 
vature in the case of the homogeneous bundle over G at 
general A 6 GL(d, R) (dropping indices a, (3), we express 
it in terms of the local coordinates A^ b by applying the 
further strain S(e x ) to I^a), and taking derivatives at 
A' = 0. This yields 

Fab.cdW = i{<PA\[Jab,Jcd]\<P\) (2.34) 
= i^ J |S(A)t[J o6 ,J od ]5(A)|vJ/>. (2.35) 

For general A, this is simply the same as at A = /, up 
to the linear transformation by A or A -1 (because of the 
unitary strain operator). This is an important general 
fact about homogeneous bundles: G/K is a homogeneous 
space (its geometry is covariant under left translation), 
and the same holds for properties of the homogeneous 
bundle over G/K (or over G). Hence we can concentrate 
on A = I. Then 

F ab , cd (0) = 5ad{Jcb) - S bc (J ad ), (2.36) 

where the expectation value is in \ipi). Specializing fur- 
ther to the G/K case with d = 2, in which ipi is an eigen- 
state of the generator of K with eigenvalue (or spin) — s, 
we have 

F ab ,cd(0) = \s(8 ad e cb - 5 bc e ad ), (2.37) 

where e ab is the two dimensional e symbol with £12 = 1. 
The information in this fourth-rank tensor can be most 
compactly expressed using the language of differential 
forms: the curvature 2-form is 

F = \ Fab,efd\ a b A dX e f (2.38) 

abef 

= - ^sidXn - d\ 22 ) A (d\ 12 + d\ 21 ) (2.39) 

at A = J. Either expression shows that if we choose 
for either ab or cd the antisymmetric combination corre- 
sponding to the generator of K, then F vanishes. So F is 
only non-zero for ab and cd in the remaining directions, 
corresponding to traceless symmetric matrices, which lie 
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along the coset space G/K . This means that the adia- 
batic connection and curvature are well-defined when we 
pass from the manifold G to G/K. 

It should be clear that the calculation is the same for 
any homogeneous bundle over G/K for G = SL(2,R), 
given a state \<fi) that is an eigenstate of the generator 
of K with eigenvalue —s. Hence it applies also to many- 
particle systems. 

6. Two space dimensions without magnetic field — paired 
states 

To obtain the Hall viscosity, we need to introduceperi- 
odic boundary conditions as in the work of Refs. (HOJIil) 
as well as vary the metric. The resulting bundles are 
still over G or G/K, but are not homogeneous bundles, 
though we will see that they can be approximated as such 
in some limits. The difference in the physical situation 
from the homogeneous bundles considered above can be 
illustrated by considering a generalization of the previous 
single-particle example. For A — I, where the metric is 
in canonical form g = /, we can take the system to be 
a (hyper-)cube (or square for d = 2) of side L, oriented 
with its sides parallel to the coordinate axes. For general 
metrics g, this becomes a rhomboid (or parallelogram for 
d = 2) if viewed in the coordinates X. The norm-square 
on the Hilbert space is now J d d x\ip(x)\ 2 where the in- 
tegral is over the d-dimensional (hyper-)cube. We will 
construct a bundle over G/K with fibre isomorphic to a 
one-dimensional space V$. Then the choice of a vector in 
each fibre can be written 

<p A (x) = /a(A t x). (2.40) 

To satisfy the boundary conditions, we require 

/a(A t (x + R))-/ a (A t x) (2.41) 

for R = Hn u n 2 ,-,n d = L(ni,n 2 , ■ ■ ■ ,n d ) with m, n d 
integers. We see clearly that /a cannot be a fixed func- 
tion /, but must have explicit dependence on A. Con- 
sequently the left action of G does not reduce to fixed 
differential operators as it did before, and so this bundle 
is not a homogeneous bundle. For the family of Hamil- 
tonians H^, the corresponding form is now 

H A = --^-V 2 x + U A (X), (2.42) 

ZlTlp 

in which ?7(X) now has explicit dependence on A, be- 
cause of the periodic boundary conditions. Sec Appendix 
[X]for some further remarks. 

We can still, however, impose an analog of the condi- 
tion that ip i be an eigenvector of the generator of K (in 
d = 2). We can require that 

^ AO (x)=e !S VA(x) (2.43) 

for O a rotation by 9. This means that /a must obey 

/ AO (0 T A T x) = e- 9 / A (A T x). (2.44) 



It is then still the case that p A for A's that differ by the 
right action of K are the same vector in the Hilbert space 
(up to a phase), so represent the same state. In eq. (|2.14|) . 
if /i corresponds to the 9 direction, then \d^f) — is\p) 
and then = 0. Consequently, there is no obstruction 
to simply identifying the fibres over points A that differ 
by right multiplication of A by O € K (up to multiplica- 
tion by a phase), and viewing the bundle as being over 
G/K instead of over G. 

Similarly, in order that f A with these properties be an 
eigenstate of H A for all A, we require that H A , and in 
particular U A , have the corresponding invariance prop- 
erty: 

C/ AO (0 T A T x) = C/ A (A T x), (2.45) 

which means that H\ is invariant under simultaneous ro- 
tation of both X and the lattice defined by the periodic 
boundary conditions. Thus the explicit dependence of 
U\ on A is in fact only a dependence on AA T , or in ef- 
fect on the metric. This deals with an issue the reader 
may have noticed, that we stated earlier that varying 
H\ with respect to A gives the stress tensor, provided A 
enters only through the metric. While the present single- 
particle model does not directly relate to Hall viscosity 
(momentum is not conserved if the potential term is non- 
zero), this question becomes relevant in the applications 
that follow, in which there is instead an interaction po- 
tential that preserves translation invariance. Varying A 
in a (translation-invariant) Hamiltonian with these prop- 
erties is the correct definition to study the effect of strain 
on the system, and ensures the symmetries of the viscos- 
ity rjabef under a •<->• b and e -f-> /. 

A particular class of functions /a obeying the condi- 
tions can be obtained by summing a function /(x) (as 
before) over translations, provided the sum converges in 
H. (We consider only d — 2, though d > 2 is similar, 
and in this case we can assume / is an eigenfunction of 
rotations.) Then we have 

/ A (A T x) = /(A T (x + R„ in2 )). (2.46) 

ni ,712 

Now if / has compact support, then for all A such that 
the support of /(A T x) does not overlap that of its trans- 
lates /(A T (x + R ni „ 2 )) for any (rii,n 2 ) ^ (0,0), the 
adiabatic curvature calculation is identical to that for 
(^a(x) = /(A T x). Thus for such /, the curvature on 
the bundle is the same as above over a portion of G/K, 
though possibly not for the more extreme A which can 
cause / to overlap with its translates. In the limit in 
which the ratio of L to the diameter of the support of 
/(x) goes to infinity, this portion of G/K becomes all of 
it. More generally, if /(x) does not have compact sup- 
port, but decays rapidly in x, then as the ratio of L to 
the decay length of / goes to infinity, the curvature in 
the bundle approaches that for the group case as ana- 
lyzed above. 

The preceding discussion can be immediately extended 
to (symmetric or antisymmetric) functions of many vari- 



ables X; that satisfy the same periodic boundary con- 
ditions in each Xj. For example, for two particles (a 
single Cooper pair) in a translation-invariant state with 
Px — iVy symmetry, the wavefunction is a function of the 
relative coordinate x = xi — X2 only. Then the pre- 
ceding discussion applies directly, as the center of mass 
coordinate drops out. It is very similar for the many- 
Cooper-pair case, viewed in real space. Because of an- 
tisymmetrization for fermion wavefunctions, this is not 
simply the sum of the single-Cooper-pair contributions 
just mentioned; there are cross-terms between different 
pairings. If the Cooper-pair "pairing" function (called 
g in Ref. Q) which appears in the real-space wavefunc- 
tion is a sum over translates of a rotationally-covariant 
pairing function (corresponding to / above) that decays 
rapidly with distance, as it does in the strong-pairing 
phase [32], then the "diagonal" (not cross-term) contri- 
butions give rapid convergence of the adiabatic curvature 
to the group-bundle result proportional to the spin — s, 
which is —1 for each pair, or —1/2 for each fermion. For 
the general terms, including cross-terms, the overlap inte- 
gral between two pairings described by two permutations 
P, P' of a reference pairing such as (12), (34), . . . , con- 
tains integrals that are a product of integrals over vari- 
ables in each cycle in P~ X P' . If these integrands have 
negligible contributions from overlaps with translates of 
the underlying rotationally-covariant function, then the 
adiabatic curvature is determined by the total spin once 
again. In Ref. Q, the treatment surrounding equations 
(2.49) to (2.51) of that paper used Fourier space, and the 
assumed form of g^ corresponds to the sum of translates 
of a rotationally-invariant g in position space via Pois- 
son summation. (The reasoning using k space resembles 
that used here in position space.) The precise expres- 
sions contain |gk| 2 /(l + |ffk| 2 ) which deals with the over- 
laps between the permutations P, P'. [The approach is 
extended a little further in Appendix 151] In the strong- 
pairing phase, convergence of the result there to that 
given by total spin only is exponentially fast because the 
expressions are analytic in k. In the weak-pairing phase, 
there is non-analyticity at small k, and convergence is 
slower, due to the long tail of the rotationally-covariant 
pairing function in that case. 

This reasoning, and the previous calculations, lead to 
the result for the viscosity tensor for a paired state in 
two dimensions (at A = 7) 

Vabcd = T^iSadCcb - S bc e a d), (2.47) 

where Stot is minus the eigenvalue of total angular mo- 
mentum in the state. Defining rf^ = flx^ix (in the ther- 
modynamic limit of a homogeneous system) we have 

r)W = ±snh, (2.48) 

where s = lim^oo s tot /N, n = ]im L ^. 00 N/L d (with 
d = 2 here) , and we have restored h to exhibit the correct 
dimensions (with s dimensionless) . The Hall viscosity is 
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minus one-half times the orbital spin density. For ^-wave 
pairs, s — —1/2; for example, for p— ip pairing, s = 1/2. 
The form (|2.48[) was obtained in Ref. [3], and holds gen- 
erally in d = 2 (though the meaning of s must be clarified 
when there is a magnetic field), and for certain compo- 
nents in d > 2 with some caveats, as we will discuss. The 
generality of the result for gapped quantum fluids, and 
the quantization of T]^ A \ will be discussed later also. 

The result for Hall viscosity was given at A = I only. 
However, as the bundle is well approximated by a homo- 
geneous bundle (under conditions that were discussed), 
the adiabatic curvature and the viscosity tensor can eas- 
ily be found at general A. We have 

^bed _ _ v (A)( g ad e bc + g bc e ad^ ^ ^ 

in which we have restored indices to their proper posi- 
tions to emphasize that this is now a covariant expres- 
sion (and is symmetric under exchange of a with b, or of 
c with d), and rf A ' is still given by eq. (|2.48|) . In par- 
ticular, while our tensors are usually written relative to 
the x a coordinates, if we use instead the X a coordinates 
in which gab = &abi 

then we see that the tensor is inde- 
pendent of the "strain" , that is the aspect ratio imposed 
by the boundary conditions. This means that the Hall 
viscosity in X space is independent of the shape of the 
system, as should be the case for a local property of a 
fluid. The form of these expressions applies to any case 
of a homogeneous bundle over G/K for d = 2. Also, 
the curvature tensor F has a similar form, by lowering 
the b and d indices using the metric g. This defines a 
G-invariant 2- form on G/K that is unique up to scalar 
multiples. 

7. Magnetic field case in two dimensions 

Next we consider particles in a uniform magnetic field 
in two dimensions. It will be worthwhile to spend some 
time on the single-particle problem. We will use the 
same gauge choice (relative to the x variables) even if 
g changes, so that we can freely take overlaps of state 
vectors even for different g. (The material in the remain- 
der of this paragraph and in the next is standard, but 
is reproduced here for the reader's convenience.) First, 
in the infinite Xx, X2 plane, it is convenient to use the 
symmetric gauge with Ax = — \Bx%, A2 — hBxx, that 
is A a = —^BeabXb, where B is the magnetic field. Then 
the covariant derivatives that act on the particle's wave- 
functions are D a = d/dx a — iA a , and [D a ,D ] = —iBe a b- 
In the conventional choice of units in which the magnetic 
length is 1, B = 1. Beginning with the metric g = I, the 
Hamiltonian for a single particle with no other potentials 
is 

H =2^2>" < 2 - 50 > 

y a 

where n a = —iD a is the kinetic momentum of the par- 
ticle and m p is its mass (the canonical momentum is 
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Pa = —id/dx a ). The kinetic momenta commute with two 
other combinations of x a and 7r Q , which can be taken to 
be the guiding center coordinates w a , 



U'a 



(2.51) 



the commutation relations with X a . The Hamiltonian in 
general is therefore 



H A = --^-Y / g ab D a D b 

2m„ 



(2.62) 



in gauge-covariant form, and have commutation relations 
[w a ,w } — —ie a b- (An alternative choice is to use the gen- 
erators of magnetic translations K a = — J2b e ao Wb, which 
obey [w a ,K ] = iS a b, [i"a>-^&] = 0-) I n the symmetric 
gauge, these are 



1 



l \^b 

O 



(2.52) 



These can be written in terms of complex coordinates 
z = x\ + iX2 as 
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Dz 


d 

dz 
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We define two sets (adjoint pairs) of simple harmonic 
oscillator raising and lowering operators by 



&t = 



,t - 



Dz' 



V2 

i 

7T 



UK 



These satisfy 



[&,&*] = [o,ot] = l 



(2.57) 
(2.58) 
(2.59) 
(2.60) 

(2.61) 



and [a, b] = [at, 6] = 0. We emphasize that the covari- 
ant derivatives, guiding center coordinates, and a and b 
operators can be constructed in any gauge, though the 
expressions in terms of x a and p a will vary. 

In terms of these operators, the Hamiltonian becomes 
H = ^-(b^b + i), and the normalized eigenstates can 
be written in terms of the normalized ground state 
cf>o = e~i\ z \ /v27r of both oscillators, atfto — bcf>o — 0, as 
(b r ) n (a r ) m ct)o/V^M.. In the lowest Landau level (LL), 
consisting of states annihilated by b, w acts in the sym- 
metric gauge as multiplication by z, and w acts as dif- 
ferentiation (times two) of the resulting polynomial in z 
that multiplies 0o [HI- 

Linear transformations in G — SL(2, R) map x to X = 

D 1 
D 2 



A T x, and so also map D = 



where g are the elements of the inverse metric to g, 
g^ 1 = A T_1 A _1 . In terms of operators on x, the gen- 
erators of these left translations are self-adjoint linear 
combinations of 



1 t 2 1 h2 
-a 1 b , 

2 2 ' 



2 2 

6+6. 



aJa - 



(2.63) 



to A : D to preserve 



The relative minus signs ensure that x a , w a transform 
by A T while n a , K a transform by A -1 . The combination 
a^a—tfb is the conventional rotation generator or angular 
momentum, and has integer eigenvalues. 

Now we can imitate the previous adiabatic transport 
calculations in the presence of a magnetic field. Working 
first in the plane, we make the (overly-naive) assumption 
that we have a single state of the form eq. (|2.32p (now 
with J's given by the expressions above; the modifica- 
tion to obtain the J generators of GL(2, R) has to be 
considered carefully as the magnetic field is not invariant 
under dilatations, but we omit details) that is an eigen- 
state of rotations of X (for example using an eigenstate 
of Ha that includes a potential term U (X) with U (X) a 
function of X 2 only) . Then the same reasoning as in zero 
magnetic field shows that we have a homogeneous bun- 
dle, and the adiabatic curvature is proportional to the 
(orbital) angular momentum eigenvalue s. [This result 
does not require that the state /(x) lie in a single LL; in- 
deed in general a potential term mixes the LLs.] This re- 
sult can be immediately generalized to any finite number 
of particles (the generators of GL(2,R) simply add); the 
Hamiltonian Hi could contain translation and rotation 
invariant interactions, and a rotation-invariant confining 
background potential term. The state is assumed to be 
a non-degenerate eigenstate of Ha- (For suitable Ha, 
this could be one of the usual trial wavefunctions in the 
fractional QH effect, in which all particles are confined 
to the lowest LL, the wavefunction is an eigenstate of 
total angular momentum, and the particles cover a disk 
in the X plane.) The resulting bundle is a homogeneous 
bundle by construction, and hence the Berry curvature 
is determined exactly by the total angular momentum, 
which is an integer. This differs from the results in Refs. 
P, 0] (which used the torus geometry) , in that the total 
angular momentum here scales as 0(i^ _1 N 2 /2) [y is the 
filling factor), which is not even extensive. We will see 
below that the difference can be traced in part to the step 
of dealing with the degeneracy of states in the plane by 
either ignoring it, or removing it with the particular form 
of potential used above. Thus the present case, with a 
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magnetic field, possesses more subtleties than the earlier 
zero-field cases. 

First, to gain insight, we return to the single-particle 
problem and examine the representation theory of the 
Lie algebra sl(2, R) implied by the above generators. We 
see that the transformations act on the a and b oscilla- 
tors separately. The single-particle Hilbert space T-L can 
be viewed as the tensor product of the two correspond- 
ing oscillator Hilbert spaces. For each such oscillator (we 
write the expressions for a only), we have the generators 
a 2 /2, a* /2, a'a+1/2, the commutators of which close on 
themselves, and so we have an infinite-dimensional repre- 
sentation of the sl(2, R) algebra. The raising operator (or 

S + ) /2 increases a^a by 2, as in the zero magnetic field 
case earlier. The representation is reducible, and splits 
into irreducible representations consisting of the states of 
even and odd a* a respectively. These are lowest weight 
representations, with the lowest S z = a^a + 1/2 values 
being 1 /2 and 3/2, respectively. If, in view of the commu- 
tation relations, one thinks of the w a as coordinates on 
phase space, then it is natural to think of the Lie algebra 
as that for the symplectic group Sp(2, R) of symplec- 
tic (or linear canonical) transformations, that preserve 
the antisymmetric form e a b- Sp(2,R) is isomorphic to 
SL(2,R), and has the same Lie algebra. However, the 
appearance of half-integer values of S z implies that the 
reducible representation constructed from the oscillator 
is only a projective representation, or can be viewed as a 
representation of the double cover SL(2,R) of Sp(2, R) 
or SL(2, R) [analogous to the finite-dimensional spin rep- 
resentations of SO (3), which are actually representations 
of the double cover SU(2)]. This representation is some- 
times called either the oscillator, the metaplectic, or the 
Segal-Shale- Weil representation [24| . 

Above we considered the bundle defined by a strain 
operator applied to a single (non-degenerate) state. To 
make contact with earlier work [l], |9[, especially that of 
Levay, we need to consider transport of degenerate sub- 
spaces. We will consider what happens when we take Vq 
to be the subspace containing all states in a given LL, and 
show that the adiabatic connection is projectively flat, 
that is the curvature is proportional to the identity oper- 
ator within the LL. This corresponds to considering the 
single-particle Hamiltonian H\ (with no potential term) , 
cq. (I2.62[) . This contains the inverse metric, and the 7r a 
or b operators, but not the a operators, and each LL is 
a degenerate subspace. The proper definition for adia- 
batic transport, and for adiabatic response, is in terms 
of varying the Hamiltonian with respect to the metric 
(see Section Hi A 1[) . and so the construction of a bundle, 
or of the action of G, is not a free choice at our disposal. 
In the present example, the Hamiltonian contains b, b^ 
only, and (for a given LL) the tensor factor in Ti of os- 
cillator states generated by a and a) is the degenerate 
subspace we wish to transport. Consequently we must 
view sl(2, R) as acting by only the b terms in the gener- 
ators, dropping the a terms, as the latter only mix the 



LL states among themselves (such operators appeared in 
Levay 0]). That is, we are free to choose a basis for the 
space V\ at each point A to be the a^a eigenstates. Then 
as varying A corresponds to Lie algebra transformations 
acting on b and b^ only, it is clear that the connection 
and the curvature are proportional to the identity ma- 
trix. [Alternatively, we mentioned earlier that if the sub- 
space Vq that is transported is the whole of T-L, then the 
curvature vanishes. The result here can be viewed in the 
same way, using the fact that the Hilbert space is a tensor 
product, and we transport all of one of the tensor factors. 
This approach shows that a choice of the action of G on 
the a, variables, or of how the orthonormal basis vec- 
tors \ip a (X)) depend on A, makes no difference in this 
case.] Moreover, the curvature must be proportional to 
the eigenvalue of (minus) the sl(2, R) generator b <! b+ 1/2 
at A = I. This gives exactly the result found by Levay 
by a similar method (on the torus), which is JV + 1/2 for 
the A/"th LL. Compared with the case of transporting a 
single state, the contribution of a* a has dropped out to 
leave this part. We will explain this result in yet another 
way when we address fractional QH states below. 

We emphasize that 6^6+1/2 is the angular momentum 
associated with the cyclotron motion on a circular orbit. 
Multiplied by u> c (uj c is the cyclotron frequency, which 
becomes l/m p in our units), it is also the Hamiltonian 
H. 

Explicit formulas can be obtained easily. If the sub- 
space Vo is the lowest LL, then a basis (unnormalized) 
is (a^) m (f)o (m — 0, 1, ...), which span the space an- 
nihilated by D/D~z. If we parametrize A g GL(2,R) 
generally by 




(2.64) 



(t 2 > 0), where 8 is a rotation angle for the element of 
K = SO(2) and I is a real number, then for the repre- 
sentative of A e G (mod K) in which 9 = I = 0, we have 

Z = (x\ + tx 2 )/t^ 2 (r = n + ir 2 ), and [l[ 

H A = ^ — [\T\ 2 Dl- Tl {D 1 D 2 + D 2 D 1 )+im^) 
2m p T 2 

= ±(-2JLJL + l). (2 . 66) 

m p V DZ DZ 2 J 

Thus for g — AA T , the corresponding lowest LL is an- 
nihilated by D/DZ, where Z = X\ + iX 2 - In terms of 
the operators b, this is solved by a Bogoliubov transfor- 
mation. Then we have a basis for the space V\ given by 
(<r) m acting on the normalized state 

(i-HY/M^Vo 

where 

a=— -. (2.68) 

i + T y ' 
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(Here \a\ < 1, which corresponds to r in the upper-half 
complex plane, r 2 > 0.) 

Passing to the case of the torus, periodic boundary 
conditions are imposed by requiring that the states be 
invariant under magnetic translations implemented by 
e _zK ' R "i "2 (or more generally, invariant up to a phase); 
this is possible only if the number of flux quanta piercing 
the square is an integer N^,, that is L 2 = 2ttN^, in our 
units. As these operators involve only the operators a, 
cr , for each LL, the boundary conditions select a finite- 
dimensional subspace of the a, oscillator space, and the 
resulting Hilbert space H of functions that satisfy these 
conditions still has the tensor product structure, in which 
one factor is finite dimensional with dimension N$, and 
the other is the oscillator space for b, W . Instead of a, a\ 
there is still an algebra of operators e - 4K R ni .n 2 / A '« that 
preserve the boundary conditions, and commute with &, 
6L Consequently, any of the arguments we used for the 
case of the plane when transporting a LL also applies 
for the same on the torus. The adiabatic connection for 
each Af is consequently the same as in the plane, and 
is projectively flat. That is, using the non-redundant 
parametrization of A by t\, T2, 0, I, and at A = /, which 
corresponds to r = i, we can change variables in the 
general result eq. (|2.39l) . and obtain for the curvature 2- 
form (removing the identity matrix in the Landau-level 
variables) 



F 



Af+1/2 



dri A dT2- 



(2.69) 



Because this example is a homogeneous bundle case, we 
can use the uniqueness of the invariant (under the action 
of G on G/K by left action) 2-form on. G/K (discussed 
above in tensor form) to deduce that anywhere in the 
upper-half plane (i.e. in global coordinates t±, t 2 ) one 
has 



F 



Af + 1/2 dn Adr 2 



(2.70) 



This was Levay's result Q , obtained by a similar operator 
point of view; there is no need for explicit reference to 
elliptic theta functions. For the A^-particle version, one 
works in the (anti-)symmetrized tensor product of single- 
particle Hilbert spaces. The transport of such a subspace 
is obtained in the same way, and is again projectively flat; 
the scalar curvature is N times the single-particle value. 
(In all these cases no approximation is needed to obtain 
a homogeneous bundle: the bundle has exactly that form 
even for finite N.) For the case of fermions filling a LL, 
N = N$, the antisymmetrized product of the LL spaces 
is one dimensional. The curvature is —N(Af+ l/2)/2, 
which agrees with Ref. Q for JV = 0, and is extensive. 
For a fluid filling the Afth LL, we have finally 



V 



(A) 



i)n. 



If instead the lowest v LLs are filled, we obtain 



(2.71) 



(2.72) 



In these examples, the orbital spin of a particle is minus 
the angular momentum of the cyclotron motion only. For 
filling the lowest v LLs (with the "real" spin of the elec- 
trons polarized), s — v/2. 

These results for non-interacting particles can be ex- 
tended easily to non-zero temperatures. Indeed, the 
derivation of adiabatic response can be extended to han- 
dle a density matrix. In equilibrium, the relevant den- 
sity matrix is the Boltzmann-Gibbs weight e~@ H where 
f3 = l/(ksT), and T is the temperature. Because this 
gives uniform weight to subspaces degenerate in energy, 
the results for the Landau-level problem are similar to 
the preceding. The Hall viscosity can be found by sim- 
ply averaging the single-particle adiabatic curvature with 
the Fermi function, times the density. It was already re- 
ported in Ref. [3J that at high temperatures this gives 
rf A > = nkjRT/(2uj c ), in agreement with the classical 
derivation [5(. 

Moreover, for non-interacting particles the standard 
linear response approach is fully tractable, and provides 
an alternative derivation. The stress tensor is the mo- 
mentum flux, 



2m„ 



TTibTTia), 



(2.73) 



where ~Ki a are the components of the kinetic momentum 
of the ith particle. The Kubo formula gives the viscos- 
ity tensor as the zero-frequency limit of the stress-stress 
response: 



T]abcd 



= lim 



-1 



i{Hi){5 ac 5 bd + 5 ad Sbc - SabScd) 

£„,(())]) 0(t)l (2.74) 



dt i 



where Hi is the many-particle Hamiltonian (i.e. kinetic 
energy), and Q(t) is the step function. (The deriva- 
tion of this formula including the "contact" term con- 
taining Hj, which is analogous to the diamagnetic term 
in a conductivity calculation, will be discussed further 
in a separate paper (26[.) The subsequent calculation 
is very similar to that for the conductivity tensor for 
non-interacting particles in a magnetic field (though in 
that case it also goes through with interactions, yielding 
Kohn's theorem). The time dependence of the stress ten- 
sor can be found explicitly; it possesses eigen-components 
(which are linear combinations of the two traceless parts 
of the symmetric tensor) that simply precess at plus or 
minus twice the cyclotron frequency. The commutator is 
then reduced to an equal-time one, which can be com- 
puted to yield the expectation value of ^2 ia 7rf a , which is 
proportional to Hj and to the b^b + 1/2 part of the an- 
gular momentum. One finds then that the final result is 
the same as above for any equilibrium system with T > 
(more details will appear in Ref. [26j). 

We emphasize that the result for the adiabatic curva- 
ture is very different from that for a disk of fluid in the 
infinite plane, treated as transporting a single state using 
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the strain operator as we did first. In that case, the fluid 
was always a disk in the X plane. By contrast, if we par- 
allel transport an initial disk at A = in the plane using 
the connection we obtained for transporting a LL, then 
apart from the LL mixing given by the function (|2.67[) . 
the fluid remains circular in the x variables (i.e. is an 
eigenstate of the sum of ala + 1/2 rotation generators), 
and so is elliptical in the X variables. In the first case, 
which had a confining potential, the fluid was rigid and 
not even strained. If we wish to find the Hall viscosity, 
this is less physical than the second case, in which there 
is a change in shape in the X plane, exactly like the shape 
of the whole system (a parallelogram in X space because 
we used the torus) for the paired states in zero magnetic 
field. In the latter system, the individual pairs do retain 
their circular form in X space, and so the total internal 
angular momentum of the pairs is obtained; this corre- 
sponds to the effect on the cyclotron variables b, b' or 
the wavefunction (|2.67[) . 

We now consider more general many-particle states in 
which the particles are strongly correlated, such as frac- 
tional QH states, that are ground states of some Hamil- 
tonian that includes interaction terms. First, we consider 
the torus, that is periodic boundary conditions. As be- 
fore, we assume that the Hamiltonian is translation and 
(in the thermodynamic limit) rotation invariant, and that 
it has a gap in its energy spectrum above the ground 
states which survives in the limit. Because of symme- 
try under magnetic translations of the center of mass 
[271 ] . all states possess an exact degeneracy of Q when 
the filling factor v = N/N^ = P/Q (with no common 
factors in P, Q). There may be further degeneracy of 
the ground states, at least in the thermodynamic limit. 
The total degeneracy (necessarily divisible by Q [13]) is 
associated with the non-trivial nature of the topological 
phase of matter. As usual we should adiabatically trans- 
port the subspace of degenerate states of interest, and 
we will do this even if there are small energy splittings 
between them in the finite size system. To simplify no- 
tation, we will sometimes ignore the degeneracy, as in 
practice in various important situations [3j it turns out 
that the adiabatic connection is projectively flat on the 
space of degenerate states. In particular, the states in 
the Q-fold degenerate space are connected by center of 
mass translations, just like the states in a single LL for 
a single particle, and we have just seen that this leads to 
a projectively flat connection. 

For the state \ip) that we consider, we have in mind es- 
pecially the ground state of some Hamiltonian restricted 
to the lowest LL, or the same in a higher LL, with the 
lower ones all filled. Clearly these forms result from weak 
interaction strengths. LL mixing is possible, but the trial 
states that serve as a starting point should usually be of 
the form stated, because otherwise the mean spin per 
particle s is unlikely to be of the form 5/2, for S a ratio- 
nal number, as we discuss below, and this suggests that 
they do not represent a topological phase. 

One might imagine that for a state in a partially filled 



LL, the Hall viscosity would be given by the same form 
as before, with s = J\f + 1/2. However, this is not 
the case: the choice of a particular subspace of ground 
states within the LL space affects the adiabatic transport. 
This choice reflects the short-range interaction Hamilto- 
nian that produces the states, which is rotation invariant 
in X space (in the large-size limit). It does not, how- 
ever, in practice (for reasonably physical states) lead to 
the restoration of the earlier result given by total an- 
gular momentum. The Hall viscosity of a large fam- 
ily of trial states in the lowest LL, given by conformal 
blocks, was calculated in Ref. Q. We will not repeat 
the earlier derivation here. The approach used in Ref. 
3] was based on the "normalization factor" argument, 
which says that if we have a normalized (orthonormal, 
in the degenerate case) section of a bundle (embedded in 
Hilbert space) that is holomorphic in its r-dependence 
except for an overall x^-independent normalization fac- 
tor, then the adiabatic connection can be found from 
that factor. The trial functions (with certain factors in- 
cluded) were argued to be normalized using screening 
properties in the 2D plasma mapping of the Laughlin 
states, and for any given more general conformal-block 
state under the hypothesis that a generalization of that 
screening holds. The necessary normalization factors (up 
to a shape-independent factor) were found by requiring 
that the short-range behavior of the interactions in the 
plasma be independent of the geometry; also it was useful 
to discretize the uniform neutralizing background in the 
plasma as a set of small point charges. The final result 
had the general form discussed above, with 

3 = i/~ 1 /2 + ^ ) (2.75) 

where v is the filling factor, and is the conformal 
weight of the field in the "statistics sector" which is 
part of the construction in the general case [T3]. The 
right-hand side can also be termed the total conformal 
weight. For trial wavefunctions more general than con- 
formal blocks, one can find s from this by the usual tech- 
niques of particle-hole inversion (for fcrmions only, and 
this is discussed further in Sec. Ill A "91 below) . flux attach- 
ment, etc. We note that once again for a large system 
the bundle over G/K is well approximated by a homo- 
geneous bundle, though the value of Stot or s may not 
be obvious from the trial wavefunction on the torus (in 
un-normalized form) , or from a Hamiltonian for which it 
is an eigenstate. 

It was found in Ref. Q that this mean orbital spin per 
particle is related to the shift. The shift S is defined for a 
system on the surface of a sphere, through the relation of 
the particle number N and number of magnetic flux 
quanta (in multiples of hc/e in conventional units, 27r in 
ours) piercing the surface for the ground state, which is 
free of defects or excitations. The relation is given by the 
form 13. 

N (j> = v- 1 N -S. (2.76) 
It was argued [ll| that the shift originates from the cou- 
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pling of the curvature of the sphere to some sort of (mean) 
orbital spin per particle, and so is given by 



S = 2s. 



(2.77) 



With s defined from the Hall viscosity, this is exactly 
what was found in Ref. [3J. 

Now we turn to the derivation of the adiabatic curva- 
ture in the plane geometry for fractional QH states, for 
an initial disk of fluid, in order to make contact with the 
approach used in this paper, in particular with angular 
momentum and homogeneous bundles. As we have men- 
tioned, it is essential to take account of the degeneracy 
of the states. We will consider here the special Hamil- 
tonians for which lowest LL ground, quasihole, and edge 
states that are zero-energy eigenstates of the Hamiltonian 
are known. These exist for the Laughlin [2g, Moore- 
Read [1^], and Read-Rezayi [3(| states, among others. 
In addition, we assume that the interaction Hamiltonian 
can be written in terms of the guiding center coordinates 
Wi, Wi only, so that it commutes with the inter-LL op- 
erators -Ki a . This enables us to separate fully the inter- 
and intra-LL contributions. 

First, we note that the general non-Abelian adiabatic 
curvature, eq. (|2.13p . for an orthonormal set \(p a (X)) of 
states depending on A can be rewritten as 

F^.ap = i [(dv,<p a \P±\d„(pp) - (d u ip a \P±\d u ipp)} , 

(2.78) 

where P±(X) is the projection operator on the subspace 
orthogonal to the "allowed" or degenerate states, 



P 



.(A) = l-£> 7 (A))fo> 7 (A)|. (2.79) 



Thus only variation of \tp) with A that takes it out of 
the degenerate subspace (for that A) contributes to the 
curvature. 

We apply this to our usual construction of a homoge- 
neous bundle, now for a degenerate set of states, 



I PA,, 



-itrX 1 J 



(2.80) 



For the degenerate subspaces that arise in the trial states 
of the QH effect, we can assume that there is a basis 
of angular momentum eigenstates. In addition, we will 
assume there is a unique state with minimum angular 
momentum in the subspace, and take this as one of the 
basis states, written as \fi.o)- This state is the trial 
"ground state" in the familiar constructions. 

For such a subspace of states, the contributions to adi- 
abatic curvature from inter-LL operators 7Tj a (or 6j, b\) 
and from intra-LL operators Wi a (or aj, a\) in the gener- 
ators J a b decouple. As only the lowest LL is involved, the 
inter-LL contribution to the Hall viscosity is the same as 
the non-interacting part discussed above. Accordingly, 
we focus on the intra-LL contribution. Instead of the 
real components J ao , it is illuminating to go to the com- 
plex components, and then the relevant non-zero part of 



the curvature, corresponding to minus the total angular 
momentum in the zero-magnetic-field cases, is (at A = I) 
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<Pi,o 



<Pi,o (2.81) 



We have set a = f3 = because our interest is in parallel- 
transporting the ground state. In this case, J^i a i lowers 
the angular momentum and must map the ground state 

t 2 

out of the degenerate subspace, while J2i a % multiplies 
the lowest LL state (in the symmetric gauge) by J2i z i/2, 
and in all the familiar cases this lies in the degenerate 
subspace. Hence the expression reduces to 



F00(0) = ( ^7,0 



T '1 
9 O 



<PI,0 



(2.82) 



This expression is manifestly non-negative. [If instead 
we consider any states a, j3 in the non-interacting prob- 
lem, then the two operators leave the state in the de- 
generate subspace (or annihilate it), and so this part of 
Fa/3 (0) cancels completely. This reproduces the result we 
discussed earlier.] Now we reverse the order of the two 
operators and obtain 



^00(0) 




<Pl,0 



(2.83) 



-N(S - 1 - v^N) 



J2 ( 



2 2 



<Pi,o 



(2.84) 



where the first term on the right-hand side of either 
line is minus the guiding-center angular momentum, 
N(Nj, + l)/2. The last term is again non- negative, while 
the first term is negative for N > j/(5 — 1), and large in 
magnitude for large N. The first term would be the full 
result if we ignored the degeneracy of the subspace, and 
transport the single (ground) state, as we mentioned ear- 
lier. In the special case of the non-interacting problem in 
which all the lowest LL single-particle angular momen- 
tum eigenstates up to are occupied by fermions (the 
v = 1 case), af annihilates the ground state, and any 
of the above expressions apply, but vanish. 

Returning to the strongly-correlated cases, it remains 
to evaluate the last term. As aj 2 appears on the right, 
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in the symmetric gauge wavefunctions it can be replaced 
by zf/2, and similarly for the adjoint acting to the left 
[25j . The expectation value we require is thus given by a 
multiple integral. This can be obtained by replacing the 
Gaussian factor in the wavefunction of the ground state 
\<Pl,o) by 

ex p(-i s E^ 2 -iEN 2 )> ( 2 - 85 ) 

differentiating the normalization integral for this state 
with respect to a and a at a = a — 0, and finally di- 
viding by the normalization factor for the unmodified 
state. In the plasma mapping for the Laughlin ground 
state [To| , the extra term in the exponent corresponds to 
perturbing the plasma by a quadrupolar harmonic po- 
tential. (This mapping is discussed further in Sec. IIIII 
below, where some justification for generalizing the map- 
ping to apply to other trial wavefunctions is also given.) 
The required second derivative is hence a quadrupolar 
susceptibility for the finite-size plasma. The plasma is 
in a screening phase, and with a perturbing potential 
will change shape so that the total electric field inside 
the region covered by the plasma is zero. If we model 
the charge (i.e. particle number) distribution as a uni- 
form charge density of v/(2ir) inside a boundary, and 
zero outside, then it is easy to calculate the response to 
the applied potential (see App. |C|). The result is simply 
is~ 1 N 2 /2; one can see that this must be so, because for 
v = 1 we can do the calculation by operator methods, 
and it is clear that for the simple form of charge distri- 
bution assumed, the result must scale as stated. (The 
result can also be extracted from Ref. [3l|, which uses a 
related approach, however that paper does not estimate 
the subleading terms which we will require.) Using this 
in general, we then obtain for the intra-LL part of the 
adiabatic curvature, 

Foo(0)=JV(S-l/2), (2.86) 

which in conjunction with the inter-LL part, which is N/2 
for the lowest LL, yields the curvature Ns, and the result 
for the Hall viscosity is again eq. (I2.48|) . as in the torus 
geometry in Ref. Q. We note that s — 1/2 is positive for 
the class of functions under discussion. 

In this argument, we made a simplifying assumption 
for the charge density. One may be concerned about 
this assumption, in particular about whether the form of 
the charge distribution near the edge affects the result, 
presumably not at the level of terms of order A^ 2 , which 
should be as stated, but at order N, and this level of accu- 
racy was required to calculate the adiabatic curvature. In 
Appendix [C] we consider the plasma arguments in more 
detail, and show that the preceding result for the cur- 
vature is correct to sufficient accuracy for our purposes, 
up to possible errors from the edge of order 0(N 1 ^ 2 ) at 
most. 



8. Higher space dimension with zero magnetic field 

In this Section we address the generalization to higher- 
dimensional many-particle systems. To simplify the dis- 
cussion, and because of the major physical applications, 
we consider only d = 3. An external magnetic field breaks 
rotation symmetry, so we set it to zero, and consider 
only paired states, as for two dimensions in Sec. Ill A 61 
above. As in that case, the paired states can be largely 
understood by generalizing the single-particle example of 
a homogeneous bundle. In the present case, we consider 
only wavefunctions /(x) that are eigenstates of rotation 
about a single axis. The functions are thus covariant 
under a subgroup K = SO(2) (which hence is compact 
but not maximal compact) of G = SL(3,R), and it is 
natural to consider a bundle over this G/K. We take 
the reference state /(x) to be covariant under rotations 
about the z axis, with eigenvalue —s as before. Then we 
consider adiabatic transport of the states I^a), viewed 
as depending on A (readers are cautioned that a corre- 
sponding Hamiltonian would here depend on A directly, 
and not only through the metric). For the homogeneous 
bundle we obtain, by similar arguments as before, the 
adiabatic curvature at A = /, 

F a b,cd(0) = -^s^2(S ad e bce + S bc e ade )n e , (2.87) 

e 

where n = (0,0,1) is a unit vector in the z-direction, 
that is along the axis of the angular momentum. The 
same result holds in local coordinates at any A, and also 
in terms of the X components, where the metric is 5 a b, 
and in these cases n e is along the angular momentum 
vector. This tensor is not symmetric under exchange 
of a with b, or of c with d, while the stress tensor is 
symmetric because of the underlying rotation symmetry. 
Passing to the paired state of the many-particle system, 
the corresponding result (in the thermodynamic limit) 
must therefore be explicitly symmetrized to obtain the 
Hall viscosity tensor, 

v abcd = -\m^2(S ad e bce + 5 bc e ade + S bd e ace + S ac e bde )n & . 

(2.88) 

The Hall viscosity response is in the plane perpendicular 
to the angular momentum vector, as was to be expected. 
There is of course also a part of the adiabatic curvature 
that is antisymmetric under exchange of a with b, or of 
c with d. This gives a Berry phase for rotations of the 
angular momentum vector, which is familiar from spin 
coherent states as mentioned earlier. The cross terms 
between symmetric and antisymmetric under the same 
exchanges vanish for symmetry reasons in our example, 
so there are no "viscomagnetic" effects in adiabatic re- 
sponse in the present states. 

Three space dimensions also brings up the topic of 
"real" spin. For a single particle Hilbert space, the spin 
enters as a finite-dimensional vector space tensored with 
the Hilbert space of functions of position. If there is 



17 



no spin-dependent term in the Hamiltonian to constrain 
the direction of the spin, then we can consider adiabatic 
transport of the degenerate (tensor factor) space of spin 
states; the latter is independent of the strain A. As in the 
case of transporting a LL, there will then be no contri- 
bution to the adiabatic curvature from the spin degrees 
of freedom. 



9. Quantization in rotationally-invariant system 

For gapped systems (topological phases) that possess 
translational and rotational invariance (in the sense that 
we are neglecting breaking of the latter by the bound- 
ary conditions, as before), we will argue here that the 
mean orbital spin per particle s (as defined in the ther- 
modynamic limit) is robust, that is it does not change un- 
der small changes in the Hamiltonian, provided no phase 
boundary is crossed. Thus it is constant throughout a 
phase. We argue further that it is actually quantized to 
rational values. 

We may compare the situation with that for the Hall 
conductivity. In a translation invariant system, the quan- 
tization and robustness of the Hall conductivity follow 
directly, and its value is given simply by the filling fac- 
tor. (For paired states in zero magnetic field, the Hall 
conductivity is zero [32j.) The situation is less simple 
for the Hall viscosity, the connection of which with rota- 
tional invariance is more subtle, particularly for the QH 
systems, as we have seen. Another approach for Hall 
conductivity that is suggestive when translation symme- 
try is broken (say, by disorder) is to average the Hall 
conductivity of a finite system on a torus over the pos- 
sible boundary conditions fa, fa (which play the role 
of A in this case). Then the integral of the curvature 
must be proportional to an integer (a Chern number), 
because the fa , fa space is compact 0, H| ■ For the Hall 
viscosity situation, the corresponding integral would be 
over the non-compact "fundamental domain" in the up- 
per half plane (for d = 2) discussed in App. [A] As this 
is not compact, no argument for robustness is evident [l| 
(moreover, the curvature, and its integral, are extensive 
in system size, and at best it would seem we might ob- 
tain quantization of Ns, not of s). Hence we must turn 
to other approaches. 

First we present a fairly simple and direct argument 
involving rotational invariance. We assume the Hamilto- 
nian conserves particle number, and so we have a ground 
state that is an eigenstate of particle number, with eigen- 
value N. In the notation of Sec. Ill A 1[ we suppose that 
(for d = 2) Ai, A2 are two coordinates on G/K, for ex- 
ample n, T2- We also suppose that the family of per- 
turbed Hamiltonian is the unperturbed Hamiltonian plus 
(in terms of X space) 

SH A = A M / d2x U » ( 2 - 89 ) 

/i=3,4,... J 

where the coefficients A M , yu. = 3, 4, . . . , can be viewed 



as further coordinates. The operators are local, and 
viewed in terms of X variables have no direct depen- 
dence on the system size or on A, except for obeying the 
boundary conditions, and we will assume they are also 
translation invariant, and would be rotationally invariant 
if not for the boundary conditions. As the perturbation 
cannot immediately close the gap in the spectrum, there 
is some neighborhood of the unperturbed Hamiltonian 
in which we may consider adiabatic transport with re- 
spect to all these coordinates. Now the "current" / M (A) 
for /1 = 1, 2 represents the (traceless part of) the stress 
tensor, integrated over space. The stress tensor of a sys- 
tem with (local) interactions is a local operator, and we 
know that its components transform like a quadrupole in 
an infinite system, due to translational and rotational in- 
variance. The expectation of the (traceless) stress tensor 
in a ground state therefore tends to zero, and presum- 
ably will do so exponentially fast in system size, due to 
the gap. This is also true for adiabatic variation of the 
perturbation coefficients A3, A4, .... Hence F^/L 2 for 
fi = 1 or 2, v = 3, 4, . . . , will go to zero as L — > 00, 
whereas F12/L 2 is of order one. This symmetry argu- 
ment is valid at all points A, using the symmetry under 
the left action of K\, which holds even in finite size with 
periodic boundary condition (for any A, it corresponds to 
rotations of X space). We can think about this in another 
way: in the thermodynamic limit, the bundle over G/K 
[with \ v , (y > 3) fixed] is well-approximated as homoge- 
neous (the local properties of a fluid should not depend 
on its shape), and G — SL(2, R) acts on it as a symmetry 
group. The curvature F Ml/ (A) for ju, = 1, 2, v > 3, can be 
viewed as a one- form on the manifold G/K (by suppress- 
ing v) , and we are saying that it becomes invariant under 
G in the limit. But as G/K is a homogeneous space, and 
in particular isotropic, any invariant vector or one-form 
field on G/K must vanish. Now it is an identity that 
(again, we neglect degeneracy of \ip) to simplify writing) 

Vh = °> ( 2 - 9 °) 

where as usual the square brackets on the indices de- 
note antisymmetrization (in terms of differential forms, 
dF = d 2 A = 0). As F ly /L 2 = F 2v /L 2 = for v > 3 
and at all Ai, A2, this implies that d^F^/L 2 = for 
/i > 3. (Strictly speaking, F lu /L 2 — > does not imply 
diFxv/L? — > 0, but we can integrate over a cube and use 
Stokes's Theorem to obtain the result.) That is, the adi- 
abatic curvature and hence also the Hall viscosity and 
s are unchanged by the perturbations, if the system is 
sufficiently large. A converse to this argument is that 
if the perturbation does not preserve rotational invari- 
ance, then one expects that in general the Hall viscosity 
will change (and the symmetry properties of the viscosity 
tensor also change, so that in general there is more than 
one independent component in rj^ e f even in d = 2). 

So far in this argument, we assumed that the particle 
number N stays constant under the perturbation. This 
is a good assumption for incompressible fluids (as in the 
QH effect), but not for compressible ones, such as paired 
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states. In particular, if in the latter we view the chemical 
potential as a parameter, then perturbing it can change 
n. In terms of states with fixed N, this can occur only if 
energy levels cross, so that N in the ground state jumps 
by some integer (most likely, an even integer) . When this 
occurs, the preceding analysis which assumes continuous 
changes in the ground state vector does not apply, so Ns 
can change. But by connecting the state at the changed 
value of N continuously with any simpler unperturbed 
state at the same N, we expect s to be the same as for 
the previous N, as long as the system remains in the 
same phase. This is the desired conclusion. 

Finally, we may consider the same question of com- 
pressible paired states for the Bogoliubov or reduced 
Hamiltonian that is quadratic in particle number, for 
which the simple BCS paired form (discussed in App. [B]) 
is exact, and with a gap function that transforms with 
a definite non-zero angular momentum under rotations, 
such as p— ip. In these, neither particle number nor an- 
gular momentum is a conserved quantum number. Hence 
we must work over G, not G/K. Under a perturbation 
of (for example) the chemical potential, the ground state 
and the expectation value of N change continuously, how- 
ever the arguments for the case of G/K no longer apply, 
and ns can change continuously. But these model states 
physically represent the same phases that can also be 
studied at fixed N as above. In both cases we have seen 
that s is determined by the angular momentum of the 
pairing, so we ascribe the continuous change to n, and 
expect that there is no change in s. 

Now we will turn to a slightly different argument 
within perturbation theory that may give more insight 
into the mechanisms for robustness. We will consider 
the effect of a perturbation in the Hamiltonian Ha for 
which the Hall viscosity is known exactly in the absence 
of the perturbation. For example, the corresponding un- 
perturbed wavefunction could be one of the paired states 
or the conformal-block QH states. The arguments we 
will give (which are for each order in perturbation the- 
ory) are somewhat schematic at this stage, and we only 
give a sketch. 

In general, the effect of a perturbation on a many- 
particle system (or quantum field theory) can be viewed 
as adding to the original ground state other states in 
which some excitations occur. These excitations can be 
factored as distinct "linked" excitations, each of which 
has to be integrated in position uniformly over the whole 
sample (due to translational invariance). (This reflects 
the linked cluster theorem.) Then we focus on a single 
such linked excitation. This object may be thought of as 
some collection of excitations (perhaps "elementary ex- 
citations" of the ground state), times an amplitude that 
depends on the separations of the excitations. We claim 
that in a system with a local Hamiltonian (both the un- 
perturbed one, and the perturbation) and a gap in the 
spectrum, these objects are local, in each order in per- 
turbation theory. That is, the amplitude decays rapidly 
as the separation of its constituent excitations becomes 



large, with the separation at which this sets in, and the 
decay rate, independent of system size as the latter goes 
to infinity. Then rotational invariance of the system in 
the X variables implies the same for this amplitude (in 
the limit) also. In this case, the Hall viscosity acquires 
contributions from each excited object which add to that 
of the unperturbed ground state. The contribution of 
the excited object can be handled similarly to the pairs 
in the paired states discussed in Section III A 61 As the 
effect of the boundary conditions drops out as the system 
size goes to infinity (due to the claimed locality of the ob- 
ject), and the object carries no net angular momentum 
in this limit, there will be no change in the adiabatic 
curvature in the limit. 

The claim of locality is clear for the linked objects of 
first order in the perturbation. These are obtained es- 
sentially by acting with the perturbation on the unper- 
turbed ground state (times an energy denominator). The 
presence of a gap should make these local. For higher or- 
ders, the constituent excitations may be created by sep- 
arate applications of the perturbation Hamiltonian, and 
so might appear not to be near one another. However, 
we believe that the result still holds, because of the gap, 
and because of the finite propagation speed of excitations 
which holds in many systems (as expressed for lattice sys- 
tems in the Lieb- Robinson bound |3J|). Indeed, because 
we consider only systems with translational and rota- 
tional symmetry, we expect that a fully rigorous proof of 
quantization of the Hall viscosity along these lines should 
be possible. Essentially, angular momentum should be 
exchanged among the particles only locally by the per- 
turbation. The total angular momentum (or the net or- 
bital spin in the QH case) in a large region should be 
unaffected. 

A third route to proving the robustness of the Hall 
viscosity (in systems with translation and rotation in- 
variance) rests on its connection with the shift S (here 
again we focus on two dimensions only), which was dis- 
cussed above, and on the quantization thereof. The re- 
lation (|2.76|) applies to quantum fluids, not only in the 
QH effect. (For examples like the paired states with no 
intensive magnetic field, v^ 1 = 0.) It can always be un- 
derstood as an orbital spin associated with each particle, 
so S = 2s, and we expect that the same s enters the 
Hall viscosity (as we have seen in examples). In the def- 
inition of the shift, N and are integers (because of 
flux quantization for the latter), at least when the par- 
ticles are bosons or fermions. For v~ x = 0, S must be 
an integer. Otherwise, if v = P/Q, where P and Q have 
no common factors, and if integer solutions (N, N^) to 
eq. (|2.76|) exist, then multiplying eq. f|2.T6[) by P we see 
that PS must be an integer. [Conversely, given integers 
P, Q, PS with P, Q coprimc, infinitely many integer 
solutions for N, N<f, can be found.] In fact, in many of 
the well-known fractional QH states, S itself is an inte- 
ger. Examples in which S is not an integer can be found 
in the par ticle-hole conjugates of the Read-Rezayi (RR) 
states [30] at general level k and with M > 1, in the 
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notation of that paper. For the RR states, 



Mk + 2 



(2.91) 



and 2s = M + 2, where k = 1, 2, . . . , and M = 0, 

1, 2, Particle-hole conjugation generally acts on 2s 

(whether it is defined via Hall viscosity or as S) as 



2s 



1 - 2vs 



and on v as v — > 1 — v. In terms of PS, we have 
PS Q - PS. 



(2.92) 



(2.93) 



Applying this to the RR states, we note that particle- 
hole conjugation applies to fermions, for which M is odd 
(30j . and one finds that for M > 1 the resulting S is not 
an integer for values k > 1. The simplest example is the 
particle-hole conjugate (at v = 3/4) of the MR state at 
filling factor v = 1/4 (k = 2, M = 3); the shift at v = 3/4 
is S = -1/3. 

That s must be a rational number, though probably 
not too surprising, was not obvious from the original defi- 
nition using adiabatic transport. [One might expect that 
arbitrary values of orbital spin for a single particle are 
allowed because any covering group of SO(2) might be 
relevant to rotations in two dimensions within quantum 
mechanics.] Incidentally, the spin-statistics relation does 
hold for s when the trial wavefunction is a conformal 
block from a unitary rational conformal field theory Q, 
that is, in these cases S is an even integer for bosons, 
an odd integer for fermions. The (orbital) spin agrees 
with that of the hole excitations, which also obey spin- 
statistics in such cases [34| (see Ref. [ill f° r a contrary 
view). However, such a spin-statistics relation does not 
have to hold for the average s when different particles in 
the ground state have different spin values. An example 
is the case of filling v LLs with fermions: a fermion in 
the Afth LL has half-odd-integer spin Af + 1/2, but the 
average gives s = v/2. 

Because 2Ps apparently must be an integer (for v^ 1 = 
0, P = 1), and given that P is fixed, s cannot vary un- 
der small perturbations. [The shift has long been viewed 
as such a "topological property" of (non-disordered) QH 
systems, even before Ref. [Hj].] It follows that the Hall 
viscosity is robust against perturbations at fixed density. 
It would be of interest to make these arguments more rig- 
orous also. Finally, we note that there is no definite N- 
Ncf, relation when disorder is present, so the shift ceases 
to have significance, due to the loss of rotational invari- 
ance on the sphere. The same will be true for the Hall 
viscosity. 



B. Numerical tests and use as diagnostic tool 

Now we turn to numerical tests. We recall that for 
parallel (adiabatic) transport of a vector around a closed 



path, in general (for a one-dimensional fibre) the vector 
changes by the phase (in the notation of Sec. Ill A 1| 



i i> A^(X)dX^ i I F^dX^dX,, 



(2.94) 



where the integral of is over a surface bounded by the 
path, and we recall that A^(X) = i((p(X)\d t _ l ip(X)) . If we 
discretize the path and the integral in small steps, so we 
have the sequence of states \tpj) = \<p(X (j))) where X(j) 
are evenly spaced along the path (J = 0, 1, . . . , M = 0), 
then we can form the product (which is manifestly gauge 
invariant) 

M-l 

n (w+ii'ft) 

Af-l 

~ J] (l + (^lv>*A M y) + 0(*Aj))j(2.95) 

3=0 

-> jfA-nW*** (2.96) 

as the size of the steps goes to zero (M — > oo). We 
will evaluate the product numerically for circular paths 
in the (ti,T2) plane, using a large number of steps. For 
comparison with the analytical result, we note that the 
relevant integral for the curvature over a disk D of radius 
po centered at (ti , T20) in the (ri, t 2 ) plane (po < T20) is 



ch~i d,T2 



2tt 



1 



V 1 - Wt 20 ) 5 



1 



(2.97) 



We will consider states with all particles in the lowest 
LL, though we know that the results can be immediately 
adapted to the case of all particles in any one higher LL. 
Using a basis of single particle states on the torus, for 
example eigenstates of e - lK i/ N <t>^ defined independently 
of A (or t), we form Slater determinants (for fermions; for 
bosons, permanents) and label an orthonormal basis of 
these by a, to obtain a basis of TV-particle states \ip a (X)}. 
A general state for N particles in the lowest LL can then 
be expanded as 



|p(A)> = £v«(A)l¥>«W>- 



(2.98) 



Normalization of |<p(A)) implies X^al^al 2 = !• The 
ground state of some Hamiltonian that acts within the 
LL takes this form, but the coefficients must be found 
for each A. We can treat the fibre of the bundle as one- 
dimensional, because in the examples we study the de- 
generate ground states have distinct quantum numbers 
and so are orthogonal for all A. Directly from the defini- 
tions, we find that the adiabatic connection is a sum of 
the non-interacting result, coming from the basis states, 
and a part from the coefficients: 



F = 



N 



- -j— 2 + i 2^ {d Tl v a .d T2 v a - d T2 v a .d Tl v a ) 



dT\/\d,T2- 

(2.99) 
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FIG. 1: (color online) s of Laughlin states for various sizes, 
showing rapid convergence with size. Both boson (v = 1/2) 
and fermion (y = 1/3) cases are shown, r = e 1 ^ 3 at the 
center of the circular path, corresponding to hexagonal geom- 
etry. The data for each case lie very close to the horizontal 
line which is the corresponding expected result. 



In effect, the non-trival second part is calculated as if 
the basis states were independent of A, and the result for 
non-interacting particles is simply added. This is very 
convenient for numerical purposes. A similar separation 
of contributions to the phase can be made in the overlaps 
for discrete steps. 

A consequence of particle-hole symmetry should be 
mentioned here (it is also mentioned in Ref. [HI): In 
eq. (|2.99p . wc know that the last term vanishes for the 
filled lowest LL. Because particle-hole symmetry can be 
defined by conjugating a wavefunction, multiplying by 
the filled LL wavefunction (in the original and additional 
coordinates) at the same flux, and integrating over the 
original coordinates, it follows easily that the last term 
reverses sign under this transformation. This can also be 
seen by rewriting the particle-hole symmetry transfor- 
mation, cq. (I2.92[) as a transformation of 2s — 1; the last 
term here corresponds to v(2s— 1). Hence it vanishes for 
the (finite-size) ground state at v — 1/2 of a Hamiltonian 
that is invariant under particle-hole symmetry. 

In Fig. [T] we show s obtained through the above pro- 
cedure by calculating the Berry phase for adiabatically 
transporting the Laughlin ground state around a circle 
in the (ti,T2) plane. For each r, the Laughlin state 
with periodic boundary conditions is generated numeri- 
cally as the zero-energy ground state of the special pseu- 
dopotential Hamiltonian |28[ on the torus. The phase 
is divided by the integral of the SL(2, R)-invariant area 
form dri A g^/tI to obtain the coefficient N~s/2. In Fig. 
[U the center of the circle is at r = e m l z (correspond- 
ing to the hexagonal symmetry case — we write r for tq 
from here on), and the radius po and number of steps 
are shown in the Figure. We verified that for such small 
radii, the result is independent of radius, and likewise 
independent of the number of steps when it is this large 
(200 steps). The circle was used to minimize effects of 
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FIG. 2: (color online) Same as Fig. [TJ but dependence on 
r at the center of the circular path is shown. Writing r = 
Irlexpi^, the horizontal axis is 9, and the corresponding |r| 
is shown for each point. The square geometry is at 6 = 90°. 
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FIG. 3: (color online) Same as Fig. [TJ but for v = 1 (boson) 
MR state for various sizes. 



finite step size relative to the local radius of curvature of 
the path, which can be severe if the path has corners (for 
example, a square). For the Laughlin states at v = 1/2 
and v — 1/3, convergence to the values predicted in Ref. 
[3| is very rapid. In Fig. [3J tests of the dependence of 
the curvature or s on the position r of the center of the 
circle are shown for several arbitrarily-chosen values of 
r, as well as r = i, the square geometry. The results are 
seen to be independent of r for moderate sizes. Thus for 
these states and for moderate sizes, s or the adiabatic 
curvature/L 2 is independent of the shape and size of the 
system, as expected for the Hall viscosity of a fluid. 

In Figs. |3] and |H results of similar calculations are 
shown for the MR state, for v = 1 (bosons) in Fig.|3l and 
for v = 1/2 (fermions) in Fig. |4] Here the ground states 
on the torus are found as the zero-energy states of the 
three-body Hamiltonian (see e.g. Ref. [29|]). The ground 
state quantum numbers are different for even and for odd 
particle numbers; for TV even, there are three sets of Q 
degenerate ground states that can be mapped to each 
other by symmetry in the hexagonal case (hereafter, all 
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FIG. 4: (color online) Same as Fig. [fl but for f = 1/2 
(fermion) MR state for various sizes. Convergence here is 
slower than previous cases. 
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FIG. 5: (color online) s and overlap-squared with the v — 1/2 
Laughlin state as a V% pseudopotential is varied; Vo = 1. Here 
N — 10. The adiabatic curvature behaves erratically very 
near the transition. 



results are for the t — e 11 ' 3 hexagonal case). For the 
Q = 2 (fermion) case, convergence is slower but does ap- 
pear to set in by around N — 14. Convergence is also 
slower for odd N than for even N. 

Next we turn to effects of perturbing the Hamiltonian. 
First, we take the Hamiltonian for which the v = 1/2 
(boson) Laughlin state is exact, namely with the pseu- 
dopotential Vo — 1, all others zero, and perturb it by 
adding a positive V% pseudopotential. A transition occurs 
at V2/V0 — 0.35. In Fig. [SI the overlap with the Laughlin 
state, and s are plotted. The overlap does not drop much 
until close to the transition point. At the same time, s 
stays close to 1, but displays large deviations as the tran- 
sition is reached or passed. It also depends strongly on 
t for V2 larger than about 0.25 (not shown). Such ef- 
fects can be attributed to increasing correlation lengths 
(which control the rate of convergence with increasing 
size) in the vicinity of the transition. 

Finally, we examine particle-hole symmetry breaking 
effects at v = 1/2 for fermions, using interactions appro- 
priate to the first excited LL, and so relevant to v = 5/2 
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FIG. 6: (color online) s of the first excited Landau level 
v = 1/2 state as a function of the strength of an ultra- 
short-range 3-body potential. For positive values s asymp- 
totically approaches the MR value. For negative values it 
passes through the anti-Pfafnan value, indicated by the lower 
dashed line, with no sign of the formation of a plateau. 



in experimental systems. We begin with the Coulomb in- 
teraction, with finite thickness (Fang-Howard parameter 
w = 2), and with V\ increased by 0.005. This Hamilto- 
nian has particle-hole symmetry, and its ground state is 
known to have high overlap with the particle-hole sym- 
metrized MR state [35|. By particle- hole symmetry, it 
has 8 = 1/2 exactly, which does not correspond to the 
value in any obvious topological phase. When perturbed 
by the three-body interaction [36|], it will eventually be- 
come the (non-particle-hole-invariant) MR state which 
has s = 3/2. It has been argued [13] that, because the 
MR phase is not particle-hole symmetric, particle-hole 
symmetry must break spontaneously at v = 5/2, and 
both MR and its particle-hole conjugate phase C(MR) 
(which has been dubbed the anti-Pfaffian) are present in 
finite-size ground states, s in the particle-hole conjugate 
phase is —1/2. Making the coefficient of the three-body 
interaction negative, one might expect the C(MR) state 
to be found, but this is somewhat naive, as there are in 
fact instabilities in this parameter region. In Fig. [Sj we 
show s for this system, for TV = 12 particles. Positive 
three-body interaction does show signs of the crossover 
to the MR value s = 3/2, while for negative values no 
saturation at —1/2 is apparent. We conclude that the 
C(MR) phase is not seen, at least for this size. This il- 
lustrates how the value of s (or Hall viscosity) can be 
used to distinguish topological phases. 

We propose that numerical measurements of the mean 
orbital spin per particle s (or adiabatic curvature, or Hall 
viscosity) can be used as a tool to determine the correct 
value of the shift. We recall that when a particular filling 
factor is studied in finite size on the sphere (or the disk 
also), one does not know a priori the value of the flux iV^ 
(resp., angular momentum) for each particle number N, 
because the shift S of the state is unknown. One may 
search for values of N, at which the angular momen- 
tum on the sphere is zero and the ground state energy 
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has a minimum or cusp (or there is a gap in the spec- 
trum). However, this is plagued with uncertainties, and 
by "aliasing", namely the fact that sequences of states 
at different filling factors can have the some N, pairs 
in common, because they have different shifts. Though 
such effects should go away as one passes to the ther- 
modynamic limit, aliasing can lead to misleading results 
given the relatively small sizes available. On the other 
hand, working on the torus provides unbiased numerics, 
as one simply seeks the ground state at N = vN^ for the 
v of interest, but in the past this gave no direct clue to the 
value of S for the corresponding state on the sphere. Our 
proposal is to find the ground state on the torus at differ- 
ent r, and evaluate s by adiabatic transport. This gives a 
numerical value for S = 2s. More importantly, the value 
of the shift is an invariant of a topological phase, so we 
gain information about the phase of matter the system 
is in. Such properties are preferable to quantities such as 
overlaps with trial states (useful though those are), which 
will ultimately tend to zero in thermodynamic limit, for 
the ground state of any local Hamiltonian except that 
which produces the trial state exactly. Here we have pre- 
sented only some demonstrations that this technique can 
work, leaving more significant applications for later work. 



III. EXACT COMPRESSIBILITY OF 2D 
SYSTEM AND STATIC STRUCTURE FACTOR 

In this section, we consider the "static" structure factor 
for ground states of QH systems, and also for some anyon 
wavefunctions. We recall that "static" actually means 
"equal time" , so that at zero temperature this structure 
factor can be computed from the ground state wavefunc- 
tion of a system. Related to the static structure factor in 
ways that we will review is a "compressibility" , not the 
physical compressibility of the 2 + 1-dimensional parti- 
cle system (which is not a purely ground-state property), 
but that of a 2D classical statistical mechanics system 
that can be defined from the ground state wavefunction 
alone. 

The motivation to consider this here came from two 
sources. The first was the version of a Hall viscosity cal- 
culation in Ref. 112], which largely follows Ref. Q, but in 
a different (cylindrical) geometry, and most interestingly 
invokes the known compressibility (in the sense just de- 
scribed) of the one-component plasma [38, 3§] related to 
the Laughlin state in order to reproduce the Hall vis- 
cosity of that state. It suggested to us that there might 
be a comparable exact result for other QH states also, 
and that these compressibilities might be (almost, but 
not quite) as robust to perturbations as we believe the 
Hall viscosity is. A second motivation was that another 
recent paper [l3| obtains a bound on the fc 4 coefficient 
of the static structure factor, that is related to the Hall 
viscosity, and it was stated that the bound seems to be 
saturated in various lowest LL trial states. 

We find exact results for this 2D compressibility, and 



hence also for the structure factor. For the QH systems, 
the result is the exact coefficient of k 4 in the small- 
wavevector expansion in powers of k. We also use the 
compressibility to recover the Hall viscosity for lowest 
LL states. 



A. Quantum Hall wavefunctions 

First we review the argument for the exact compress- 
ibility in the one-component plasma. We recall that the 
normalization integral for the Laughlin state | with 
exponent Q (filling factor 1/Q) is the partition function 
of a one-component plasma with a uniform neutralizing 
background (l0| : 



(3.1) 



j ' ~\\d 2 Zi exp 



Q^Tmlz.-^f-i^N 2 



To go further, we must make the system explicitly neutral 
by restricting the background charge to a disk of radius 
R containing charge N (in units where the particles carry 
charge 1), and also include the self-interaction of this disk 
of charge. Then we have instead 



Z = J Y\_d 2 Zi expQ 



^ln|zj - zj\ 2 
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+ dVp(x')ln|x 4 -xf 

i 

+ d 2 xd 2 x' p(x)p(x') In |x - xf 
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where p(x) = —N/(irR 2 ) for |x| < R, otherwise. This 
change should make no difference to correlations well in- 
side the boundary, providing the system is in the screen- 
ing phase, p = —l/(2nQ) inside the disk produces the 
same density as the standard Laughlin state. In the stan- 
dard notation for the one-component plasma, T = 2Q 
here. 

We are interested in the isothermal compressibility of 
this system, which is related to the response at fixed tem- 
perature (and fixed Q) of the density of particles to the 
total potential, including that produced by the other par- 
ticles, as well as the background potential. This can be 
found by a scaling argument, apparently due to Ref. [Hj]. 
The area covered by the particles is determined by the 
background charge density or the potential it produces, 
assuming that the system is in the screening phase. The 
dependence of the free energy on the area of the system 
can be found by scaling: if the area f2 = ttR 2 is changed 
to ttR 2 X 2 (note that this changes p also), then by rescal- 
ing lengths r = r A we have 

Z(X) = \2N+QN(N-l)-2QN*+QNi Z{1) (33) 
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X (2-Q)N 2il) _ 



(3.4) 



Thus the dependence of Z on S7 has been determined. 
With Z = e~ F ', and the pressure p = —{dF/d£l)p,N (we 
set the inverse effective temperature to one as it is purely 
conventional), we obtain the equation of state 



(3.5) 



where n = N/Q is the number density. The isothermal 



compressibility is then defined as xt 
and is given by [39[ 



n 1 {dn/dp) T ,N 



(3.6) 



A negative value of this compressibility does not imply 
an instability in this system, due to the long-range inter- 
action. For (5 = 0, Xtq —n is the ideal-gas result. 

This argument generalizes immediately to other states 
in the lowest LL, as follows. We see that the important 
points were: (i) the treatment of the background charge 
with uniform density that is set to —i>/(2t:) after taking 
derivatives; (ii) the wavefunction (other than the back- 
ground charge parts) was homogeneous of total degree 
NN^/2, where N^, is given by the relation (|2.76[) . For 
the Laughlin state, we had v~ x = S = Q. Then using 
the general form for a general lowest LL state we obtain 
similarly 



Xt 



V 



(3.7) 



(y drops out). This is the first main result of this Section. 
We can replace S/2 by the mean spin per particle s if we 
wish. 

We should emphasize the conditions under which this 
argument is meaningful. The result gives the isother- 
mal compressibility, which is supposed to be an intensive 
thermodynamic property of the system, provided that 
the boundary effects are negligible. This holds if the 2D 
system is in a screening phase for charge, and also screens 
(has exponentially decaying correlations) for the other, 
"non-charge" or "statistics" sector also. This is the same 
hypothesis under which results were obtained for states 
constructed from conformal blocks in Ref. Q. Note also 
that the 2D Coulomb interaction of the one-component 
plasma has been supplemented by interactions that are 
neither two-body, nor simple to write as a Hamiltonian 
(logarithm of the Boltzmann weight). We assume that, in 
spite of their long-range appearance, they do not produce 
any net long-range two-body number-number interaction 
and so do not require any neutralizing background. The 
filling factor is therefore determined by the charge sector 
only (assuming screening) ; this is a standard line of rea- 
soning in QH systems. The argument here is so general 
that it still applies directly even for a ground state that 
is obtained as a perturbation of a trial state, as long as 



it remains in the lowest LL and in the screening phase. 
Thus it is robust within a topological phase within the 
lowest LL, as long as rotational symmetry holds. 

By combining this result for the c omp ressibility with 
the derivation of Tokatly and Vignale [12| , we can recover 
the result for the Hall viscosity. In their derivation, the 
quantity that enters is the "interaction" part of the bulk 
modulus, which corresponds to x^ 1 — n in terms of the 
above. We see that, when multiplied by — h/2, this is 
precisely the Hall viscosity result, eq. (|2.48[) . This estab- 
lishes a connection between rf A > and the compressibility, 
and reaffirms the connection Q with the shift S, for LLL 
states. By comparing the various arguments, we can see 
that, apart from the different geometry, the approach in 
Ref. [12j essentially takes a different route to the same re- 
sult as Ref. Q for the normalization of the ground state, 
or at least for a relevant derivative of that normaliza- 
tion with respect to r, which is also connected with the 
derivation of xt above. 

We now turn to the "static" (equal time) structure 
factor of the system. It can be defined in terms of the 
(number) density-density correlation function at equal 
time USES: 



5(x,x') = m6n(x)5n(x')\y)/(nZ) 



(3.8) 



where Sn(x) = n(x) — (^t'ln(x) \^)/Z (note that the 
factor 1/Z is required because |\f f ) introduced above 
was not normalized). After taking the thermodynamic 
limit for fixed x and x', we can assume that S(x, x') 
is translationally and rotationally invariant, and write 
S(x — x') = 5(x, x'). 5(x) is related to the two-particle 
reduced density matrix g(x) of the state |\E r ) by 



5(x) = 5(x) + nh(x), 



(3.9) 



and h(x) = g(x) — 1 — > as x — > oo. As these defini- 
tions involve only the coordinates of the particles, and 
not the momenta which would require differentiation of 
the wavefunction they can be defined in exactly the 
same way for any classical system by replacing \^>\ 2 / Z 
by the probability density for the particles' coordinates. 
s(k) is now defined by taking the Fourier transform of 
S(x). Because we assume rotation invariance, we will 
write it as s(k). 

The 2D Coulomb interaction in the charge sector, as 
the only effective long-range force, requires the back- 
ground charge, and can be separated out in the same 
way as in the one-component plasma (see Ref. [40( for a 
detailed discussion). Then in the screening phase, stan- 
dard arguments lead to 



k 2 

*(*) = pr 



k A X^ 
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-+o(fc 4 ) 



(3.10) 



(the standard notation f(k) — o(fc 4 ) means /(fc)/fc 4 — > 
as k — > 0). Here k 2 D is the inverse Debye length squared, 
and the appearance of xt expresses the "compressibility 
sum rule" [4pJ (similarly, the vanishing coefficient of k°, 
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FIG. 7: (color online) The static structure factor s(k) of the 
v = 1 MR state for bosons, with N = 100 particles. Mea- 
surements were taken over 2 x 10 s MC steps. The inset shows 
s(k) over the entire range of k. 



and the fixed coefficient of fc 2 , are due to "charge neu- 
trality" and "perfect screening" sum rules respectively). 
For the present case, fcf, = 2 in our units (HI, and we ob- 
tain for a trial state in the lowest LL the static structure 
factor, exact through 0(fc 4 ), 



1 



1 



s(fc) = -r + -(5/2 - l)fc 4 + o(fc 4 ). 



(3.11) 



This is the second main result of this Section. For the 
Laughlin state 5 = Q, it was obtained in the same way in 
Ref. [4l| . It does not seem to have been known previously 
that the fc 4 coefficient is robust within a phase, even for 
the Laughlin case. 

We have given the result here for the full static struc- 
ture factor as this seems to us more natural. The lowest- 



LL projected structure factor is given by [41 1 
So (fc) = s(fc)-(l-e- fe2 / 2 ). 



(3.12) 



Thus for srj(fc) one should drop fc 2 entirely, and add 1/8 
to the coefficient of fc 4 , which becomes (5—1) /8. We note 
that (5 — l)/2 is s with the non-interacting or inter-LL 
contribution 1/2 subtracted off. Haldane [13| obtained 
such a relation as an inequality by analytical arguments, 
and also found that it appears to be an equality in nu- 
merical calculations for some examples, all for rotation 
invariant systems. We find that it is exact for all low- 
est LL states under conditions that should correspond 
to their being in a topological phase, as long as rotation 
invariance holds. 

We have tested the above prediction for the exact coef- 
ficient of fc 4 in the case of some trial states, using Monte 
Carlo (MC) techniques on the torus. For both the Laugh- 
lin and MR trial states, MC is tractable. We focus here 
on the MR state at v = 1 for bosons. The wavefunctions 
for these states on the torus have been found previously 
(see e.g. Refs. [29|,|32|])- We consider any one of the three 
states for N even; these are related by symmetry for the 
case of the hexagonal system r = e l7r / 3 . First in Fig. [7J 
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FIG. 8: (color online) The small k behavior of s (k) (the 
lowest-LL-projected structure factor) of the v = 1 MR state 
for 100 bosons, obtained from s(k). Also shown is the ex- 
pected small k behavior fc 4 /8. 
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FIG. 9: (color online) The coefficients of k 4 in s(k) - k 2 /2 
for the two smallest non-zero k values, plotted versus for 
various sizes. The expected value as N — > oo is shown as a 
horizontal line 



we show s(k) in full and at low fc, for 100 particles. In 
Fig. [8j we make the subtraction on the same data set to 
obtain so(fc), and compare with the expected behavior 
fc /8 (no fitting parameter). The agreement is good. 

In Fig. [9l we show results from MC at different particle 
numbers. For each size, we measured s(k) at the smallest 
two non-zero fc values only, subtracted fc 2 /2, and divided 
by fc 4 . There is a very clear trend towards the expected 
value 5/2 — 1 = in this example. 

In the preceding results for the compressibility of the 
2D system, the correction to the ideal gas law p — n (for 
P = 1) contains s, and one might expect this always to 
be the same as the Hall viscosity. However, this is not 
so. For a wavefunction with all particles in the Afth LL, 
the wavefunction contains Yii times a corresponding 
lowest LL function, plus terms in which z^s and Zj's can- 
cel in pairs (for each i). We emphasize that here it is 
not correct to interpret each d/dzi, as we are not 

in the lowest LL. The latter would be valid if we were 
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calculating the total spin of the wavefunction, as for the 
Hall viscosity. But for the compressibility of the 2D sys- 
tem, we instead need the scaling dimension, which is the 
degree under the rescaling of the polynomial part of the 
wavefunction under z — > \z, z — > Xz. The wavefunctions 
with any particles in the J\f > LLs are not in general 
homogeneous under such a rescaling. One might proceed 
by using some average degree, but its value is difficult 
to predict. Hence the results for the compressibility and 
the k A term in the structure factor are not valid in higher 
LLs, nor for LL mixing with an initial lowest LL state, 
unlike the results for Hall viscosity. 
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B. Anyon wavefunctions 

We may make a similar argument for trial states whose 
wavefunctions are conformal blocks without a neutraliz- 
ing background (e.g. by "removing the charge sector" 
from the QH functions). These can describe states of 
some kind of anyons (not necessarily Abelian). An ex- 
ample would be the SU(2) level one conformal block for 
spin- 1/2 primary fields, 



n«0fc - w i) 



1/2 



(3.13) 



which is a two-component state (zi, Wi, i = 1, . . . , A/2, 
are the coordinates of the spin f and i particles, respec- 
tively) [42j (other examples, including that corresponding 
to the p+ip superfluid, were also discussed in Ref. [3[). 
Such a function can be put in a finite area by transferring 
it to the torus, for example. The same scaling argument 
now gives Xt 1 = (•"■ — h)n, where h is the conformal 
weight of the field representing a particle in ^ (h = 1/4 
in the example), and replaces S/2 in the above (indeed, 
one may think of these wavefunctions as having v = oo, 
and S = 2h). Then for a system without a long-range 
Coulomb force or neutralizing background, one has for 
the structure factor at k — (actually defined as the 
limit k — > after A -> oo) 



ns(0) = djl 



(3.14) 



by the fluctuation-dissipation theorem applied to the 
classical system for T — 1 (/i is the chemical potential), 
and dn/dn = tT 2 xt by a thermodynamic argument. This 
yields 



8(0) 
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We should mention also that the expected Hall viscos- 
ity of the system is rf A ^ — \hnh~ [3], though this value 
might be corrected due to the presence of gapless excita- 
tions. We note that for sufficiently small h (h < 1/2 in 
a paired example like that above; in general the condi- 
tion depends on the operator product structure implied 



FIG. 10: (color online) s(k) at the smallest non-zero k for 
systems of up to 150 anyons in the wavefunction in the text. 
The horizontal line is the expected value s(0) = 4/3. 



by the wavefunction), the integral of |^| 2 is convergent 
at short separations, but for large values of h is not. In 
the latter short distance cut-off will be required, 

and it is not clear if this defeats the scaling argument. 

For stability, dn/dfi should be positive, and hence (ig- 
noring the concerns at the end of the previous paragraph 
here) h < 1. The latter is the condition that the per- 
turbation of the vacuum by the field be a relevant per- 
turbation (as discussed in Ref. Q), and in the present 
situation this is necessary in order for the system to be in 
the "screening" phase in the statistics (not charge) sec- 
tor. Thus a transition in behavior should certainly occur 
if h passes 1, where the system does not screen. (The 
structure of the wavefunction obtained from a conformal 
field theory will depend usually on h, so it generally can- 
not be varied continuously; nonetheless the point holds.) 
If h > 1 one expects the particles to form bound clus- 
ters and drop out of the long distance behavior instead 
of exhibiting screening; the structure of the clusters may 
involve a short-distance cut-off scale. For the example 
above, |\I/| 2 can be viewed as a 2D plasma with both + 
and — "charges," which correspond to the spins t and 
4,; for this plasma, the conventional parameter T = Ah. 
Screening holds for the spin density correlations, because 
in the example h — 1/4 < 1 [42]. Physically as a wave- 
function, it is a (paired) charge superfluid, and a spin 
liquid with a Hall conductivity for the spin current. The 
preceding arguments imply that it has s(0) = 4/3. 

More generally, if the wavefunction is scale-covariant 
as well as translationally invariant and rotationally co- 
variant, but not holomorphic (as is a conformal block, 
away from the diagonals Zi = Zj), the same arguments 
hold with 2h replaced by x, the scaling dimension of the 
fields in l^l 2 , and x < 2 for screening and stability. 

To test the prediction on the example, we performed 
Monte Carlo simulations on the torus, that is with peri- 
odic boundary conditions, in order to calculate s(k). For 
this we need the wavefunction on the torus, which can 
be obtained from the Coulomb plasma on the torus as in 
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Ref. [|| . The wavefunction is not single valued, instead it 
changes when one particle (of either spin) performs a cir- 
cuit around one of the cycles. This means that there are 
in fact two "conformal blocks". The Boltzmann weight 
for the plasma can be obtained, following for example 
Refs. EES!; it is 



2 |* B («i,...;ioi > ...|r)| 2 , 

e=0,l 

where the two conformal blocks are 

* e (z i ,W7 fe |r) = 
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where Z = ^J2i z i ~ {~I2k w ki z ij = %i ~ z h an d 
Wij = Wi — Wj. $ a .b{z\T) are elliptic theta functions with 
characteristics, defined by 

O)6 (z|t) = J2 e™ T{n+a)2+2m{n+a){z+b) . (3.18) 

n 

The function rj (r) is the Dedekind function (q — e 27TlT ), 

oo 

V(r) = q 1/24 l[(l-q n )- (3-19) 

n=l 

The function £ (z/L x \t), 

(3.20) 
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[where we have used the more familiar notation $i (z |r) = 
— "&i/2,i/2{ z \ T )]i ' IS called the prime form for the torus, 
and is periodic up to phase factors. Using L 2 T2 = L 2 , 
which is held fixed as r varies, we can write the conformal 
blocks as 

% e (Zi, Wk\r) = 

e/2 . o (2Z/L x \2r) n 1/2 
»7(r) f j 1, v{t) 

n / ^iK/^|t) \ 1/2 n f Mh-w k ] /L x \T) \ - 1/2 



which implies that if a Zj (wfe) is increased by L x t 
{—L x t), the two blocks are mapped to each other (by 
a unitary transformation), as was stated above. 

We have tested the prediction on the trial state by 
Monte Carlo, using these wavefunctions on the torus. We 
estimate s(0) simply as the value of s(k) at the smallest 
non-zero k. The results shown in Fig. [10] are in excellent 
agreement with the predicted value s(0) = 4/3. 



IV. CONCLUSION 

To summarize, the main results are that the general 
connection of Hall viscosity to mean orbital spin, uncov- 
ered in Ref. Q, can be understood from the purely ge- 
ometrical relation that the commutator of distinct shear 
operations is a rotation. This brings in the angular mo- 
mentum eigenvalue (in the simplest cases). The robust- 
ness of the Hall viscosity to perturbations of the Hamil- 
tonian, as long as one remains within a topological phase, 
has been shown, making use of rotational invariance. The 
relation to the shift [3] is known from all examples, and 
has been shown for all lowest LL states, but still lacks a 
truly general derivation. When proved directly, this re- 
sult will be a theorem that relates two physical properties 
defined separately. In addition, we performed numerical 
tests of the stated results, which are very convincing at 
least for trial states. Moreover, we gave a proof that 
the 2D compressibility and coefficient of fc 4 in the static 
structure factor are also related to the shift [l3j], and 
again tested these numerically on trial states, this time 
by Monte Carlo calculations with up to 100 particles. 

The orbital spin is revealed here as a true emergent 
property: it was not evident microscopically that there 
is a conserved quantity of this type. Macroscopically (i.e. 
on scales larger than the correlation length, and at en- 
ergies below the gap), there is a well-defined spin per 
particle, and hence an orbital spin density. When the 
particles are in motion, there must also be a conserved 
(spatial) orbital spin current density. 

Many open problems remain for Hall viscosity. Cal- 
culations can be extended to non-rotationally invariant 
systems. Further tests of its robustness to perturbations 
can be made. It can be used as a diagnostic for which 
topological phase a system is in. Most importantly, it 
will be interesting to find techniques to experimentally 
measure it in any of the systems discussed here. 



x [(im-r^/Vr)]^ 2 , 



up to z-, w-, and r-independent factors. This form is 
more convenient for exhibiting the symmetry properties 
(note that the Gaussian factors cancelled). We note the 
property of the theta functions in the center of mass (Z) 
factor, 

d a . {z + t\2t) = e-^ T ^"^ a+1/2fi (z\2T), (3.22) 



(3.21) 
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Appendix A: Automorphic properties 



Here we address briefly some further points in the case 
of systems with periodic boundary conditions, using the 
example from Section III A 61 though similar considera- 
tions apply mutatis mutandis in all other cases. If the 
elements of A' are integers, then A' T R ni „ 2 is again a lat- 
tice point R n ' n ' for all n\, n<i. Such matrices form the 
discrete subgroup V = SL(2, Z), which again acts on the 
left on G, as A —> A'A with A' e r. In general, this does 
not map <pa to itself, but only maps it to another state 
= ^A'A with the same properties (note that the left 
action of T commutes with the right action of K) . It is 
similar in the case of the Hamiltonians H\ : while for all 
A, L^a(X) has the symmetry under all translations of X 
by A T R ni „ 2 , it may not obey invariance under the left 
action by T, t/ A / A (X) = E/a(X) for A' e V. However, in 
many physical situations one would require such invari- 
ance, so that the system does not depend on the choice 
of the basis vectors for the lattice. Examples may be 
obtained by summing a rotationally-invariant U (X) over 
the lattice translations, similarly to the wavefunctions 
in Section III A 61 (One important special case is a peri- 
odic S- function.) Then in many cases the non-degenerate 
eigenstate /a(X) must be mapped to itself (as a func- 
tion of X) by elements of T, that is for any A' g V, 
/a'a(X) = /a(X) up to a phase and some action of K 
(this assumes that energy levels of H\ do not cross as A is 
varied to implement the transformation in T, as would be 
required in order to use the adiabatic theorem). A func- 
tion /a(X) that is invariant under T (up to a phase and 
K action) would be called an automorphic function (for 
the group of translations and operations in T); then (if / 
is also an eigenstate of K) we have a line bundle over the 
space Y\G/K, In the case d = 2, r is called the modular 
group, and T\G/K can be represented by a fundamental 
domain in the upper half plane G/K . (In cases with a set 
of degenerate functions _/a(X), one must unfortunately 
speak of "vector- valued automorphic functions".) Physi- 
cally, we should view /a'a(X) and /a(X) as representing 
the same state. However, it is important to realize that 
even when /a(X) is invariant (or invariant up to a phase, 
and so on) under the left-action of T, this does not mean 
that the state y A (x) is; instead, ^a'a(x) = ^a(A' t x). 
We can recover <^a(x) by changing variables A' T x — » x. 
Such a change of variable induces a unitary mapping of 
the Hilbert space (of functions of x) to itself if the bound- 
ary conditions are unchanged as they are in our example 
(more generally, a change in boundary condition, and a 
corresponding change of Hilbert space may be involved 
01 9)- For a non-trivial closed path in T\G/K, the latter 
transformation is a form of monodromy [3[ ; the adiabatic 
transport around the path produces a well-defined holon- 
omy (up to a phase associated with the path dependence 
due to the curvature), which gives a "holonomy repre- 
sentation" of the modular group T, which was studied in 
Ref. Q, and will not be considered further in this paper. 



Appendix B: Paired states in k space 

In this Appendix, we give an expression for the adi- 
abatic curvature for a paired state, working in k space 
with (anti-) periodic boundary conditions. Expressions 
of this type were given in 0] , but in terms of = Uk / 
(these quantities will be defined below). However, this 
form may raise concerns in the weak coupling region in 
which it might be that «k is zero over a range of k, say 
|k| < ko. In order to be certain that this does not change 
the claimed result, which leads to a quantized s in the 
thermodynamic limit, we will obtain here an expression 
that is more general and deals with this situation. 

The BCS ground state for a translationally-invariant 
spinless or spin-polarized system on the torus is 

i^=n'( uk +^ci ct -k)i°)- (bi) 

k 

Here itk and t>k are complex functions of k, and the 
product Yl is over each distinct pair (k, — k). To be 
slightly more general than elsewhere in this paper, we 
can take the boundary conditions in the two-dimensional 
system to be either periodic (+ or m — 0) or antipcri- 
odic (— or m = 1) in each of the two directions. The 
allowed values of k are then given in complex form by 

1/2 

k x +ik y — 2iTi(n2 ~ uit) / (Lt 2 ), where we parametrize 
A by t = T\ + iT2 as in Sec. Ill A 71 These values of k re- 
fer to X space, but note that in x space the plane waves 
arc independent of r or A. m and ni arc defined by 
ri\ = ra\j1 (modulo integers), and similarly for n^. For 
H — h or mi = riJ2 = 0, the value k = is omitted from 
the sum, but an additional fermion occupies that mode 
if the system is in the weak pairing phase [32| . We may 
consider any of these boundary conditions. The state is 
normalized provided |uk| 2 + |vk| 2 = 1 for all k (we may 
define itk and i>k for all k, i.e. both k and — k, in some 
suitable way). 

We may immediately find the Berry or adiabatic con- 
nection, here for complex r as the generalized coordinate, 
A T =i(ip\d T tp) (d T = d/dr), 

A T = J ^2 ("k^ U k - (<9 r Uk")w k + ^r^k - (<9 T Uk~)u k ) ■ 

(B2) 

(The sum is over the allowed k; k = is omitted in the 
++ case.) Now we proceed similarly as in Ref. 0. We 
will consider the p— ip paired state. We may make the 
gauge choice, for example, that — it(|k|) is real and a 
function of |k| only, while 

1/2 

Uk = mk|), (B3) 

n 2 - nir 

where w is real. Then 

4 ni — n\T 
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This is exactly what was found in Ref. [3[ (with the same 
gauge choice), but here we have a more general deriva- 
tion. From this point, one can differentiate with respect 
to r (to obtain the curvature) either before or after tak- 
ing the thermodynamic limit in which the discrete k sum 
becomes an integral. The resulting adiabatic curvature 
leads to precisely the same Hall viscosity or s @ as dis- 
cussed in the main text, from a different point of view. 
This is a consequence of the assumed form for Uk, i>k 
(though independent of the choice of gauge), which are 
rotationally-covariant functions of k, and depend on r 
only through the discreteness of k. [The present trial 
state is not an eigenstate of particle number N, unlike 
those in the text, but standard arguments about the size 
of number fluctuations in it suggest that if we project 
the state to definite N, we obtain the same result with 
N replaced by N.] The finite-size result may be differ- 
ent for other assumed forms, or for the state resulting 
from solving the gap equation in the finite geometry. We 
have not investigated these possible effects further, but 
we note that the general quantization argument in Sec. 
Ill A 91 suggests that, in a physical setting, the result for 
s in the thermodynamic limit of a rotationally-invariant 
system cannot change unless a phase boundary is crossed. 



Appendix C: Quadrupolar susceptibility of a 2D 
plasma 

Here we consider in more detail the quadrupolar sus- 
ceptibility of the 2D system, as discussed in Sections 
III A 71 and IIIIl For the Laughlin state, the 2D system 
has partition function as in eq. (|3.f P (in this Appendix 
we use the uniform background charge density of infinite 
extent, the potential of which is shown in this equation). 
The modification of interest is to add to the exponent 
the harmonic potential acting on all the particles, 

-5«il>i-v?), (ci) 

i 

where we took a = a = u\ real for definiteness. For trial 
states other than the Laughlin state, there are additional 
interactions in the plasma, but we assume these do not 
change the essentials of the plasma arguments, as they 
only change short-range correlations (see Sec. . 

To estimate the free energy of the plasma in the pres- 
ence of the perturbation, or its second derivative, the 
susceptibility, we should use the correlations of the par- 
ticle density. With the long-range 2D Coulomb inter- 
action between the particles, we expect that the main 
effects can be handled as a "self-consistent" field or po- 
tential produced using the mean number density as the 
source. Technically, this means that in the correlation 
functions, we express them in terms of 1-interaction irre- 
ducible parts, and the latter are short range (this is like 
the random-phase approximation, but here we cannot as- 
sume translational invariance, because of the boundary 



of the plasma) . [This is also part of the standard argu- 
ments used in the bulk of the system in Sec. IIIIl [40].] 
Little or no information about these short-range effects 
will be needed, other than the assumption that screening 
holds. 

We will treat the plasma from a macroscopic view- 
point, using electrostatics, but in order to deal with pos- 
sible microscopic effects and length scales (related to the 
short-range correlation effects), we will later include pos- 
sible densities of higher multipoles, in addition to the 
charge density. In the screening phase, the plasma be- 
haves similarly to a conductor. In standard treatments 
of electrostatics [45| , a conductor is viewed as a charge- 
neutral system with fixed boundaries, but with mobile 
charges in the interior. In an external field, a charge 
density accumulates on the surface (edge, in 2D) so that 
the net electric field in the interior vanishes. The sur- 
face charge per unit length of edge is proportional to the 
discontinuity in the normal electric field. By contrast, 
in the plasma we have no a priori position for the edge, 
which can move in response to the field. Hence the (self- 
consistent) edge is not only an equipotential, so that the 
tangential field vanishes, but also the normal component 
of electric field must vanish, as does the surface charge 
density. We note that, because our perturbation is har- 
monic, the interior of the perturbed plasma has the same 
constant charge or number density as the unperturbed 
one. At this level of treatment, the charge density is 
viewed as constant in the interior, with (in the macro- 
scopic viewpoint) a step at the edge. 

To find the shape of the drop in the presence of the 
harmonic perturbation, we can consider the (real and 
single- valued) potential experienced by an additional test 
particle outside 51, which can be expanded in the form 

oo 

$(x) - ^\n\z\ 2 -\\z\ 2 +Y.^rnZ- m + ® m ~ m ), (C2) 

m— 1 

where <E> m are constants. The coefficients can be related 
to the multipole moments of n(x), in a standard way 
[4o| . The charge density inside Q does not change, so the 
potential there is given by a similar expansion (without 
the first two terms, and with positive, rather than neg- 
ative, powers of z and z). The net change in potential 
due to the perturbation must be independent of x inside 
f2, and this condition determines the change in shape f2. 
In polar coordinates (r, 9) , we find that the radius as a 
function of polar angle 9 is 

r(6) = R(1- a x cos 29), (C3) 

to first order in ai; here R is the unperturbed radius, 
irnR 2 = vR 2 /2 = N. Thus for a small quadrupolar har- 
monic perturbation, Q has a small elliptical deformation. 
Then the potential outside is 

$(x) = NQ\n\z\ 2 -\\z\ 2 -\a l {x 2 -y 2 )+\a 1 ^{x 2 -y 2 ) 

(C4) 
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to first order in a%, and up to an additive constant. The 
curve on which — V$ = for this <E> reproduces eq. (|C3[) . 

In reality, the edge of the QH trial state is not sharp, 
but rounded on a scale of order the particle spacing (or 
a screening length in the plasma). We can include this 
in the macroscopic viewpoint to a first approximation as 
a (normally-oriented) dipole layer on the surface, due to 
moving some charge in- or outward, compared with the 
step distribution. In the presence of a dipole layer, the 
potential has a discontinuity on crossing the layer, equal 
to the dipole moment per unit length D. However, for 
a conductor, or for our plasma, the potential both just 
inside and just outside the edge must be constant along 
the edge, and hence the dipole moment is independent of 
position on the edge (though its magnitude may depend 
on the shape, in principle). The charge density of our 
plasma can then be modeled to this approximation as 

n(x) = ^-I n + DV 2 I n + ..., (C5) 

Z7T 

where In is the indicator function for the region in 
which the charge density is non-zero, that is Iq = 1 in- 
side, = outside. The existence of such a dipole layer 
can be seen by evaluating the expectation of £^ \ z i\ 2 m 
the unperturbed plasma, and also by viewing it as related 
to the angular momentum in the unperturbed wavefunc- 
tion. For this state, we evaluate: 

E^W 2 ) = l N ( N + + 2 ) ( C6 ) 

= -v~ l N 2 + 2irv~ l DN + ...,(C7) 

where the first equality comes from the operator view- 
point and the value of the angular momentum, and for 
the second we used eq. (|C5I) to evaluate the left-hand 
side. The terms of order TV 2 agree, while the terms of 
order N give 

O-^l, (OS, 



for the circular edge. This seems to be related to the 
dipole moment on the edge discussed by Haldane jl3j . 

Then the quadrupole moment of the first term in n(x) 
in eq. (|C5|I should still agree with that in the potential. 
This is the case, because the second term gives zero: 



D / d 2 x(x 2 - y 2 )V 2 I n = (C9) 



by integrating by parts. (Similarly, the contribution of 
D to all other multipole moments is zero.) Hence D 
plays no role in the quadrupolar susceptibility, which is 
determined only by the shape of f2, to our accuracy. For 
consistency, the dipole moment should not contribute to 
the electrostatic energy either. This is the case because 
it enters as the integral of DE n along the edge, where 
E n is the normal electric field, and E n is zero, even for 
the deformed shape. Consequently, we see no reason why 
D should change in the presence of the perturbation, as 
the local environment at the edge is unchanged. There 
may be effects involving more derivatives which do cause 
D to change weakly with the perturbation, however it 
does not contribute to the quadrupole moment in any 
case. Higher multipole moments at the edge involve more 
derivatives and do not contribute either. We conclude 
that the susceptibility required in the main text is simply 



-v~ l N 2 +0{N 1 ' 2 ), (CIO) 



The error term is present because there may be true edge 
effects in which the free energy of the plasma depends on 
the length of the edge due to short-range effects that we 
have not considered. Under the small elliptic deforma- 
tion, the change in the perimeter of VL is of order the 
unperturbed perimeter, times a 2 , which gives the possi- 
ble error term of the order as stated. 
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